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Abstract 

We consider instanton solutions of Euclidean Hof ava-Lifshitz gravity in four dimen- 
sions satisfying the detailed balance condition. They are described by geometric flows 
in three dimensions driven by certain combinations of the Cotton and Ricci tensors 
as well as the cosmological-constant term. The deformation curvature terms can have 
competing behavior leading to a variety of fixed points. The instantons interpolate 
between any two fixed points, which are vacua of topologically massive gravity with 
A > 0, and their action is finite. Special emphasis is placed on configurations with 
SU(2) isometry associated with homogeneous but generally non-isotropic Bianchi IX 
model geometries. In this case, the combined Ricci-Cotton flow reduces to an au- 
tonomous system of ordinary differential equations whose properties are studied in 
detail for different couplings. The occurrence and stability of isotropic and anisotropic 
fixed points are investigated analytically and some exact solutions are obtained. The 
corresponding instantons are classified and they are all globally R x and complete 
spaces. Generalizations to higher-dimensional gravities are also briefly discussed. 
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1 Introduction 

Based on ideas that were originally developed in condensed matter physics lUl and 
later applied to the description of aspects of particle interactions lEHlUlHSl, a modifi- 
cation of general relativity was recently proposed in IIU. |Z1 and further studied under 
the name of Horava-Lifshitz gravity. In this theory, which includes higher-order cur- 
vature terms on spatial slices, the diffeomorphism invariance of general relativity is 
broken explicitly setting a privileged time direction. This affects the ultra-violet be- 
havior of the quantum theory, which, hence, looks power-counting renormalizable. 
Some efforts have been made to prove consistency of the quantum theory ||8l l9l, but 
the number of propagating degrees of freedom seems to invalidate the matching with 
Einstein's gravity in the infrared regime, and hence seems to disprove this theory as a 
viable alternative to general relativity l|T0l[m[T2l[T3| . The investigation of these issues 
is still going on. Despite the difficulties and reservations, the Hofava-Lifshitz gravity 
still provides an interesting classical and quantum field theory framework, where one 
can address some interesting questions and explore several connections to ordinary 
gravity or string theory. These also include the appearance and relevance of geomet- 
ric flows, which is the main subject of the present work. 

Geometric flows, and, in particular, Ricci flows are interesting in their own right. 
In mathematics they turned out to play a crucial role in implementing Hamilton's pro- 
gram for proving Poincare's and Thurston's conjectures | |T4l[T5| (but see also fTE] and 
| [T7| and references therein). In physics they originally appeared in off-critical string 
theory via the renormalization-group equations of two-dimensional non-linear sigma 
models, where the evolution of the metric under the Ricci flow equations provides the 
running of the bulk coupling to lowest order in perturbation theory (see IITSlI for the 
original result). In this context, the renormalization-group time is provided by the log- 
arithmic length scale of the world-sheet, but in some cases it can also assume the role 
of genuine time, describing real-time evolution in string theory in regimes where the 
friction due to the motion of the dilaton effectively reduces the second-order evolution 
equations to the first-order renormalization-group flow equations ||T9||20|| . 

Ricci flow models also appear in the framework of four-dimensional gravitational 
instantons of general relativity. Solving Einstein's equations is, in general, an impos- 
sible task. It is substantially simplified under the assumption of self-duality as a suf- 
ficient condition to find vacuum solutions in the Euclidean sector of the theory. Ho- 
mogeneity of spatial sections is often a further simplification to find explicit solutions. 
Although what we call space is somewhat arbitrary in Euclidean gravity, the latter 



2 



statement can be made precise by assuming a foliation in three-dimensional leaves 
that are invariant under an isometry group of motions. For these particular vacuum 
solutions, it turns out that the Euclidean time evolution of the homogeneous leaves 
inside the gravitational instanton can be recast as Ricci flow equations for the corre- 
sponding geometry on the homogeneous model spaces ||2ni20ll22| . 

The modification of gravity proposed by Horava in ||6l[7|l shares some features with 
the previous setting that allow to define the analogue of gravitational-instanton con- 
figurations. In particular, a foliation of the four-dimensional space is assumed from 
the very beginning with a privileged time direction at the level of the action. Fur- 
thermore, a condition called detailed balance, which is borrowed from non-equilibrium 
thermodynamics, requires that the dynamics follows from an appropriately chosen 
three-dimensional superpotential action. In the Euclidean version of the theory, this 
resembles the self-duality condition with similar consequences: for a class of configu- 
rations that minimize the action, the time evolution becomes first-order and describes 
a geometric flow on the leaves of the foliation. Unlike general relativity, where Ricci 
flow is equivalent to self-duality only for configurations with homogeneous leaves, 
the description of instanton-like solutions by geometric flows is generic in Hofava- 
Lifshitz gravity. The nature of the corresponding flow depends on the choice of the 
three-dimensional action used for detailed balance and it is driven, in general, by a 
certain combination of Ricci and Cotton tensors as well as the cosmological constant 
term. It is not our concern, and it will not be pursued here at all, to find whether 
such combinations of curvature tensors can also arise from the renormalization- group 
equations of some quantum field theory. 

Our aim in this paper is to investigate aspects of the Euclidean dynamics of Hof ava- 
Lifshitz gravity with detailed balance. Our motivations are diverse. First, classical 
instanton-like solutions are important for the determination of transition amplitudes 
in quantum gravity. They are also useful in the Hartle-Hawking formulation of quan- 
tum cosmology 12311 , even though classical cosmology per se requires the analysis of 
real-time equations. The geometric flows that emerge in this framework involve ten- 
sors with higher-order spatial derivatives terms, and, as such, they are new in the 
literature; they reduce to previously studied examples only for some special values 
of their parameters. Thus, it is instructive to formulate the flow equations in all gen- 
erality, determine the nature of the fixed points and their stability properties, obtain 
explicit solutions, as well as study general questions such as the monotonicity of the 
evolution, the possible formation of singularities, the occurrence of bounces and so on. 
These questions arise naturally in the general theory of geometric flows and they are 
bound to be relevant for the space-time interpretation of the analogue of gravitational 
instantons in Hof ava-Lifshitz theory. 
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The answer to these questions will be accomplished partially using some ansatz 
for the underlying three-dimensional spaces, leading to mini-superspace truncation 
of the flow equations. Otherwise, it does not be seem possible to draw general con- 
clusions for the general system of equations, at least at the current level of our un- 
derstanding of this problem. Still, the results that will be described are indicative 
of what should be expected in general. More systematic investigation of the infinite- 
dimensional dynamical system at hand requires substantial mathematical work that is 
not contained in this paper. Following the paradigm of gravitational instantons with 
isometry groups in Euclidean Einstein gravity, we will consider homogeneous geome- 
tries on the three-dimensional spatial slices of Horava-Lifshitz gravity and focus, in 
particular, to the case of Bianchi IX geometry as a class of homogeneous but generally 
non-isotropic deformation of S'^; this model is often referred to as mixmaster universe 
in the Lorentzian (real-time) approach to cosmology ||24ll25ll26l , and was recently dis- 
cussed in the framework of Hofava-Lifshitz gravity ||27ll28| . Other Bianchi classes as 
well as more general inhomogeneous deformations of the three-sphere (under appro- 
priate ansatz) can also be studied along similar lines, but they will not be discussed. 

We set up the general problem using the Bianchi IX model geometry and study in 
detail some specific examples of the flow for different couplings in Hofava-Lifshitz 
gravity. Even in this case the resulting equations in mini-superspace are not easily 
tractable for generic values of the couplings. First, we will consider the Ricci flow 
and some of its variants that describe solutions of the modified theory of gravity with 
anisotropy scaling parameter z = 2 (see next section for this and other definitions) 
and compare them to instantons with SU{2) isometry in ordinary gravity. We will 
also consider the Cotton flow, separately, and use it to construct solutions of Horava- 
Lifshitz theory with anisotropy scaling parameter z = 3 by ignoring all Ricci curvature 
terms that become subdomtnant when the volume of is very small. We will also 
consider the combined Ricci-Cotton flow and explore the equations in detail first in 
the limit that the speed of light vanishes or equivalently Newton's constant becomes 
infinite (it is often called Carroll limit after [|29|) . The normalized Ricci-Cotton flow 
and the unnormalized variant of it with vanishing cosmological constant provide the 
relevant equations in this limit. Finally, we will consider the general Ricci-Cotton flow 
with arbitrary couplings and obtain several qualitative results for its solutions. 

In all cases it is assumed that the parameter of the superspace metric of the theory 
is restricted to values A < 1/3 (in which case the cosmological constant will also be 
taken non-negative) so that the flow equations extremize the classical action, up to 
important boundary terms. Proper account of the boundary terms leads to the defi- 
nition of instanton solutions as finite-action trajectories that interpolate between fixed 
points. Our analysis shows that the Ricci and the Cotton tensor terms can compete 
with each other, and, depending on the relative sign between the two, the flow equa- 
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tions can exhibit symmetric as well as anisotropic fixed points. The nature of these 
fixed points and their stability properties also have implications for the space-time 
interpretation of the corresponding gravitational instanton solutions. Axisymmetric 
solutions are associated with spaces with SU{2) x U{1) isometries, and, hence they 
are easier to describe in closed form. 

Instanton solutions of Horava-Lifshitz gravity, as they are defined, are rather spe- 
cial configurations that rely on the existence of multiple-degenerate vacua and cor- 
respond to special flow lines, which guarantee finiteness of their action. As it will 
turn out, they are also free of singularities and their space-time metrics are regular 
and complete. Note, however, that other flow lines, possibly with infinite action, also 
describe solutions of the second-order equations of motion, but they may have singu- 
larities. Although we are primarily interested in the class of instanton solutions, one 
should be open-minded for other more general possibilities too. For this reason, as 
well as for mathematical completeness, we will investigate the phase portraits of the 
flow equations in all generality. The selection of special trajectories that correspond to 
instantons will be made much later together with their space-time interpretation. 

In section 2, we first briefly review the formulation of Hof ava-Lif shitz gravity with 
detailed balance condition putting emphasis on the structure of the potential term and 
its associated superpotential. This analysis is then carried to the Euclidean regime 
where "zero-energy" (i.e., self-dual-like) configurations exist satisfying the flow equa- 
tions. Restrictions on the parameter A are also obtained together with an entropy func- 
tional that changes monotonically along the flow lines. The results are then used to 
define instanton solutions as in ordinary point particle systems. The Bianchi IX model 
geometry is introduced in section 3 where the truncation of the Ricci and Cotton flows 
are studied separately in detail. Section 4 is entirely devoted to the analysis of the nor- 
malized Ricci-Cotton flow and the explicit construction of its axisymmetric solutions. 
Section 5 discusses the case of unnormalized Ricci-Cotton flow obtained for general 
couplings. It contains the case study of positive and zero cosmological constant for 
axisymmetric configurations with A < 1/3. Section 6 is devoted to the space-time 
interpretation of the flow line, as gravitational instantons, making comparisons with 
the analogous instanton solutions arising in general relativity. Complete classification 
of all gravitational-instanton metrics with SU(2) isometry is also obtained. In section 
7 we outline generalizations of the framework to higher-dimensional Horava-Lifshitz 
gravities, which, for instance, in 4 + 1 dimensions give rise to a new system of flow 
equations on four-manifolds driven by the Bach tensor. Finally, section 8 contains our 
conclusions and poses several questions for future work. 
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2 Non-relativistic gravity, detailed balance and flows 



2.1 Non-relativistic gravity: a reminder 

The theory of non-relativistic gravity developed in ||6HZl is valid for general space- 
time dimension D + 1. It has three main features: 

• Space- time is assumed to be topologically Aio+i = IR x A^o, leading to a nat- 
ural codimension-one foliation. Diffeormorphism invariance is broken down 
to the subgroup of foliation-preserving transformations. This breaking is con- 
trolled by a parameter A. 

• Scaling properties of space and time are different and captured by an integer z. 
Power counting renormalizability of the theory requires z = D. 

• The interactions are determined by a detailed balance condition following from 
a Euclidean D-dimensional diffeomorphism-invariant action, which gives rise 
to marginal and relevant terms in D + 1 dimensions. 

The last item above is not generic in Horava-Lifshitz gravity and it can be relaxed 
by allowing more arbitrary coefficients for the various marginal and relevant terms. 
However, it is a necessary ingredient in our study to connect it naturally with the 
theory of geometric flows. Thus, detailed balance will be assumed in the following. 

Let us adopt the ADM (Arnowitt-Deser-Misner) decomposition of the metric (see, 
for instance, ||30| '), which is suitable for the D + 1 foliation of space-time, 

ds^ = -N^dt^ + g,j [dx' + N'di) [dxj + NMi) , (2.1) 

where N' and N are the shift and lapse functions respectively. Here, i, j,. . . run in D 
dimensions and all tensors that will appear in the following are D-dimensional. 

Using this decomposition, the Einstein-Hilbert action in D + 1 dimensions (up to 
total derivative terms that may contribute in topologically non-trivial spaces) reads as 
follows, 

Seh = J d°+ix v^N (K.jK'j _ ^2 + _ 2a) , (2.2) 

where A is the genuine cosmological constant in D + 1 dimensions. In this expression, 
Kij is the second fundamental form that measures the extrinsic curvature of the leaves 
at constant t, 

^^i = ^ i^^Srj - V,N^- - VjN,) , (2.3) 



6 



and its trace K = g^^Kjj is the mean curvature. The first two terms in equation (2.2 1 
provide the kinetic energy, since they include time derivatives of the field gij. Their 
specific combination can be recast in the form 

^G'S^^u = K^jK^' - (2.4) 
using the DeWitt metric in superspace 

G^=l(gV + gV]-g^^g''. (2.5) 



DW= 2 [S +S g' ) -g 'g 

The potential term in Einstein gravity is provided by the three-dimensional Ricci 
scalar curvature R and the four-dimensional cosmological constant term A (when it is 
present), as shown in ( Z2] >. 



In non-relativistic gravity, space and time scale as [t] = — z, [x] = —1 and it is 
further assume c^that [Ni] = z - 1, [N] = and [gij] = so that [K^] = Iz; it should be 
contrasted to general relativity where space and time scale the same with z = 1. This 
asymmetry is further implemented in the action, both in the kinetic and the potential 
terms by requiring foliation preserving covariance. The kinetic term is generalized as 

Sk = ^ j dt dPx ^NKijG'^^^Ku = ^ / d^x ^gN (^KijK'j - \K^^ , (2.6) 

where [k^] = z — D and A is a dimensionless coupling measuring the breaking of the 
full diffeomorphism group. Here, G'-'*^^ is the generalized metric in superspace 

G'^"' = \ [fg^' + g''^^) - Hh^'- (2-7) 

that coincides with the DeWitt metric when A = 1. It is worth stressing that this metric 
can be positive-definite or indefinite depending on A. Indeed, gu is an "eigenvector", 

G'^'^'gu = (1 - . (2.8) 

with eigenvalue 1 — AD. The sign of the latter changes at A = 1/D where the inverse 
no longer exists. Thus, for A < 1/D the metric is positive-definite and it becomes 
indefinite for all A > 1/D that include, in particular, A = 1. This behavior and the fact 
that A is ultimately an unprotected parameter of the theory, which, in principle can 
take any real value, should be kept in mind when considering quantum corrections. 



^Proper restoration of the speed of light, which scales as [c] = z — 1, explains the various dimensions, 
as described in detail in the original works |l6l[7|. 
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The potential term of the theory has the general form 

Sy = -|dfd°x^N%] (2.9) 

and can also contribute in various ways to the breaking of diffeomorphism invariance. 
Note that [V] = z + D and there is a large freedom to choose V so that it includes oper- 
ators of dimension less than or equal to z + D (called relevant and marginal operators, 
respectively). In order to reduce this freedom and take advantage of the renormal- 
ization properties of a D-dimensional system, it was proposed in [6l 13 to introduce a 
detailed balance condition that allows to express the potential in terms of a "superpo- 
tential" as follows: 

V = "^E^^GrjuE"', (2.10) 



where 

and W a D-dimensional action so that [E'^] = D. The tensor Qiju is defined as 



(2.11) 



1 A 

QijU = 2 {Sikgjt + gUgjk) - J^j^^gijgM (2.12) 

and coincides with the inverse of the metric in superspace with generic A, i.e., 

G'j'^'Gurnn = li^U + . (2.13) 

The resulting theory is not invariant under general coordinate transformations of 
space-time. Indeed, since A4d+i is topologically R x A^p, it is only appropriate to 
consider invariance of the action under the restricted class of foliation-preserving dif- 
feomorphisms, 

t = t{t), f = x\t,x). (2.14) 

Then, the lapse function N associated with the freedom of time reparametrization is 
restricted to be a function of t alone, whereas the shift functions N, associated with dif- 
feomorphisms of A^d can depend on all space-time coordinates. This is often called 
the projectable case of Horava-Lif shitz gravity and it will be assumed in the following. 
The non-projectable version of the theory leads to dynamical inconsistencies Il311 . 

The choice of W depends on the dimension D. Here, we recall the choice for D = 3 
with z = 3 that ensures power-counting renormalizability of the four-dimensional the- 
ory; generalization to higher dimensions will be discussed later in section 7. Then, the 
marginal operators in question are obtained from the three-dimensional gravitational 
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Chem-Simons action, which is familiar from topologically massive gravity ||32ll , 

Wcs = — / cv3{cv) , (2.15) 



with density given in terms of the connection one-form (X) by 



073(0;) = ^Tr l^o; Ada;+ ^o; Ao; Ao; ) . (2.16) 



The corresponding variation gives 



where C^^ is the Cotton tensor and e^^^ = 1. The Cotton tensor is traceless, conserved 
and it vanishes identically for conformally flat metrics. 

Relevant operators in four dimensions are generated by the Einstein-Hilbert three- 



dimensional action 

Wi 



EH = ^ / d°xyg(R - 2Aw). (2.18) 

Note that neither is the four-dimensional Newton's constant nor Aw is the four- 
dimensional cosmological constant, but they will be identified shortly. The variation 
of this action leads to 



Combining the Chem-Simons and Einstein-Hilbert contributions to E^^ , with their 
respective couplings, the full potential of Horava-Lifshitz gravity reads 

2 2 



V = -^C^Qj ^C'R^j + 



f^ur. 4A-1 j\ 



w 

.2 A 

W 



The ultra-violet behavior of the resulting theory is dictated by the quadratic Cotton 
curvature term, which is marginal with dimension 2z = 6, and corresponds to z = 3. 
It improves a lot the ultra-violet behavior of ordinary Einstein gravity at the expense 
of breaking general covariance of the theory at short distances. When the Cotton term 
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is absent, the resulting theory has a potential with quadratic Ricci curvature terms that 
become dominant in the ultra-violet regime and so z = 2. In either case, in the infrared 
limit one expects to flow by the most relevant operators (of dimension 2 and zero), 
which correspond to the last terms in equation ( 2.20[ >, and recover general relativity 



provided that A also flows to 1. However, no rigorous proof of any of these statements 
is yet available in the literature. Also, the counting of physical degrees of freedom of 
the theory, which is crucial for viewing it as viable modification of general relativity, 
is obscured by the outcome of local invariances and their potential restoration. 

We also recall for completeness, using the relativistic coordinate = ct, that the 
effective speed of light for general A is given by 



k2 / A 



^ (2.21) 



2k2,V1-3A 

with [c] = z — 1. This shows that Aw must be negative when A > 1/3 to ensure 
reality of c; likewise. Aw must be positive when A < 1/3. Also, the four-dimensional 
effective cosmological constant is given by 

A = ^Aw (2.22) 

and, therefore, the range A > 1/3 does not allow for de Sitter-like backgrounds in 
Hofava-Lifshitz gravity. These identifications are necessary in order to compare the 
infrared limit of the deformed theory to ordinary gravity so that the effective Newton 
constant reads as ^ 

Gn = ^ . (2.23) 
32nc 

Furthermore, A should approach (flow to) 1 in the infrared limit in order to recover 
the full reparametrization invariance of general relativity, but this particular problem 
will not concern us at all here. 



The search for classical solutions requires the use of the potential (2.20) and it is 
impossible to solve in full generality. Symmetry ansatz such as spatial homogeneity 
makes the problem more tractable, but still not exactly solvable. This includes, for 
example, the case of Bianchi IX geometry leading to the mixmaster universe model in 
four space-time dimensions with Lorentzian signature. We will not pursue this line of 
investigation here (see | [27| for a detailed analysis and comparison with the mixmas- 
ter universe in general relativity ||24l |25l ), but elaborate on the Euclidean version of 
Hofava-Lifshitz gravity and then analyze its instanton solutions for Bianchi IX spatial 
geometries. 
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2.2 Euclidean action and flow equations 



Besides the various physical motivations pertaining to the analysis of the Euclidean 
version of Horava-Lifshitz theory, there is also a technical advantage for constructing 
solutions that satisfy first-order equations in time. This possibility is also encountered 
in general relativity when self-duality on the Riemann (and more generally on the 
Weyl) tensor is imposed leading to gravitational-instanton solutions in the Euclidean 
regim^ Although there is no direct analogue of self-duality in gravitational theories 
with anisotropic scaling, the detailed balance condition offers the appropriate replace- 
ment for defining instanton-like configurations. This is in fact possible in all dimen- 
sions unlike gravitational instantons of ordinary gravity that are only defined in four 
space-time dimensions. 

The Euclidean action is obtained by setting t — > —it, iW , whereas iS is 

traded for — S. Using equations ( Z6] > and ( Z9] > one obtains 



j dt dPx v^N (^^ [K.jK'j - AK^) + . (2.24) 



The expression for the potential (2.101 allows to rewrite the Euclidean action (2.241 in 
the form S = S' + S" , where 

S' = ^Jdt d°x v^N (^K,j ± y e^^-^^^nE-"^ G^j'^' (^K^e ± y 6^fc£.sE") (2.25) 

and S" is a total-derivative contribution to the action ||6l. This boundary term will be 



considered later in detail (see equation (2.331). The different signs correspond to the 
choice of time direction. 



The action (2.251 is bounded below by zero provided that the superspace metric 
Qtjke positive-definite. Then, configurations that obey the first-order differential 
equations 

Krj = ^\QrjU^^', (2.26) 

are extrema of the action and as such they provide solutions of the Euclidean theory S'; 
the leaves of the corresponding space-time foliations have prescribed extrinsic curva- 
ture. This possibility arises only when A < 1 / D, in which case we must also demand 
that A^v is non-negative so that the speed of light in the Lorentzian version of the the- 
ory is real. Otherwise, for A > 1/D, the action is non-bounded below by zero and 



^Self -duality is best described in terms of the curvature two-form, as TZah = ^T^ab ^ orthonormal 
frame ||33ll . These equations are second-order in time, but they can be integrated once to yield first- 
order equations [34 1 that will be paralleled to the rnstanton solutions of Hofava-Lifshitz gravity; for 
further details see also [22 1, where this analogy is made even sharper for gravitational instantons with 
homogeneous spatial sections. 
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the first-order equations are not guaranteed to provide classical solutions. Thus, from 
now on, we restrict ourselves to A < 1 / D and > 0, where the first-order equa- 



tions (2.26) provide extrema of S'. They are also extrema of the action S provided that 
the boundary term S" is properly accounted. This problem will be treated carefully in 
section 2.3 and lead to the precise definition of instantons. 

The solutions that we will investigate can be expressed in the form of geometric 
gradient flow equations for the metric gij modulo reparametrizations generated by 
the shift functions, 

^tg^j = ^K^NGiju^^^ + ^t^j + V;N, . (2.27) 

Since we are only considering the projectable case of Horava-Lifshitz gravity, N is 
only a function of t and can be absorbed by redefining time, as N{t)dt — dt. It is 
also natural to define vector fields with components = Ni/N that generally depend 
on space and time coordinates. Then, the geometric-flow equations assume the more 
standard form that will be used in the following. 

Specializing to D = 3, we write down explicitly the flow equations obtained by 
combining the variation of the Chern-Simons and Etnstein-Hilbert actions. 



l-^ f 2A-1 Ayj \ k"^ ^ „ ^ „ ^ , ^ 



choosing for definiteness one of the two sign options; the other follows by time rever- 
sal. These equations describe the parametric evolution of a family of three-dimensional 
geometries 

ds^ = gij{x;t)dx'dx^ (2.30) 

driven by the Ricci and Cotton tensors and the cosmological constant term, and, as 
such, they will be called Ricci-Cotton flow equation^ 

The fixed points are determined (modulo reparametrization terms) by the solu- 
tions of three-dimensional topologically massive gravity: 

-^Q/ = R„- - 2Awg„- , with R = 6 Aw • (2.31) 

^cs 



■'Perhaps a more appropriate name is Ricci- Yamabe-Cotton flow, since Rgij is the driving term of the 
so called Yamabe flow, dtgij = —Rgij. The latter was introduced in the literature [35J to solve Yamabe's 
conjecture stating that any metric is conformally equivalent to a metric of constant curvature. Its effect 
is complementary to the Cotton term of the flow, which changes the conformal class of the metric. Thus, 
the combined flow equations we have obtained contain several competing deformations of the metric. 
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They include Einstein metrics with vanishing Cotton tensor, like the round sphere 
metric on (for Aw > 0), which is homogeneous and isotropic. There are other fixed 
points, however, with constant scalar curvature but with non-vanishing Cotton tensor. 
We will see later, as example, that they correspond to particular homogeneous but 
non-isotropic metrics on S^. The coexistence of fixed points from different conformal 
classes of the metric make this flow particularly complex. 

The driving terms of the Ricci-Cotton flow involve, in general, third-order deriva- 
tives in space coordinates (originating from the Cotton tensor), and, therefore, it is not 
possible to apply standard results from the mathematics literature to prove even the 
short-time existence of solutions. Nevertheless, the mini-superspace models that will 
be studied later show that these flow equations are well-behaved and the trajectories 
converge to fixed points after sufficiently long time. 

Some special cases are worth noting, since they have already appeared in the liter- 
ature for different reasons: 

• ^cs ~^ the Cotton tensor contribution drops out and one obtains a variant of 
the Ricci flow on three-manifolds, which is second-order and well studied in the 
literature. Its trajectories describe solutions of z = 2 Hofava-Lifshitz gravity in 
3 + 1 dimensions, whereas the fixed points are Einstein metrics -R^ = 2Aw gij. 

• — ^ oo: the Ricci and cosmological constant terms drop out and one obtains 
the pure Cotton flow that was recently introduced in the literature ||36|. Its trajec- 



tories describe solutions of z = 3 Hofava-Lifshitz gravity in 3 + 1 dimensions, 
in the limit under consideration. The fixed points are conformally flat metrics, 

Crj = 0. 

Even these simpler cases are impossible to solve in all generality. Mini-superspace 
models have been used to study the long time behavior of the Ricci and Cotton flows 
for homogeneous geometries H37ll36B . 

Apart from the Ricci and Cotton flows that will studied separately in the next sec- 
tion, there are also some other special cases of the combined Ricci-Cotton flow that 
are relatively easier to study. First, by considering that limiting case A — > — oo, which 
lies in the allowed range A < 1/3, one obtains the normalized Ricci-Cotton flow, 
which is driven by a traceless tensor, and, thus, preserves the volume of space; the 
cosmological constant decouples in this case. It becomes relevant in the Carroll limit 
of Horava-Lifshitz gravity, where the effective speed of light vanishe^ The effective 

'^When A ^ — 00, the metric in superspace becomes singular as its inverse has zero eigenvalues. Yet 
the flow equations are well-defined and so is the potential term Sy of the gravity action. It is opposite 
to the case A = 1 /3 for which the metric in superspace has zero eigenvalues and its inverse becomes 
singular; the latter case corresponds to the limit of infinite speed of light, where the theory develops 
anisotropic Weyl rnvariance. 
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speed of light vanishes also when Aw = 0, irrespective of A, and the corresponding 
flow equations will be studied separately. The general case, with arbitrary coefficients, 
is much more complex. The pattern of fixed points and specific trajectories will only 
be discussed for axially symmetric deformations of S^, which correspond to solutions 
with SU{2) X U{1) isometry 



2.3 Entropy functional and action bound 



From now on, and in all examples that will be studied later, we consider flows without 
the effect of space reparametrizations, setting Nj = 0. Also, we will take advantage 
of time reparametrizations to set N{t) = 1 for convenience. We will also assume that 
the spatial slices are compact spaces without boundary. Here, we provide an entropy 
functional for the geometric flows arising in Horava-Lifshitz gravity in arbitrary di- 
mensions. This functional is also be related to the lower bound of the Euclidean action 
S (rather than S') when boundary terms S" are properly taken into account. 

When the metric in superspace is positive-definite (choosing A < 1 / D in D spatial 
dimensions), the superpotential functional W changes monotonically along the flow. 
This follows easily by considering 



dW 

"dT 



(2.32) 



which is the integral of a quadratic quantity, and, therefore, increases or decreases 
monotonically depending on the overall sign. Thus, W is an entropy functional for 
the flows under consideration. 

With this in mind, let us revisit the original Euclidean action S of the theory and 
its lower bound taking into proper account the boundary terms. Equations ( 2.24[ > and 
( |2l0l > yield 



S = ^ I dtd^'x^K.jG'i'^Kke + y | dtd^'x^WQi^uE 



di d°x ( ^ ± ^^.^.^E-") a^^' (ku ± y 



T2 / dtd^'x^gK.jE'K 



(2.33) 



The first term in (2.33 ) is S' , given in equation (2.25 1, and the last term is the advertised 
boundary contribution S" . For positive-definite superspace metric, the Euclidean ac- 
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tion S is bounded from below by the boundary term S", because S' > 0. Thus, 

S > t2| dtd^x^K.jE^j = T j dtdPx^E'jdtg.j = ±\ j dt^, (2.34) 



having set N(f) = 1. The time integral of equation (2.32 1 shows that the lower bound 
of S is always positive, as expected. 



The flow equations (2.26) provide time-dependent extrema of the action S' . They 
are actually its ground states, since they make S' vanish. Since S" is a boundary term, 
these ground states are also extrema of S under appropriate boundary conditions that 
make the variational problem well-posed. This can be easily verified for the class of 
flows with finite action (i.e. finite S"). Note for this purpose that the fixed points of the 
flow are static solutions of both S and S' , since they are, by construction, critical points 
of the D-dimensional action functional W sitting at the minima of the Horava-Lif shitz 
potential. If different minima exist, they will be all degenerate with zero potential en- 
ergy. Hence, time-dependent solutions that interpolate between any two fixed points 
are guaranteed to satisfy the equations of motion following from the Horava-Lifshitz 
action S. These solutions have finite action, given by the value of the boundary term, 
and is natural to call them instantons as they interpolate between two different static 
minima, which are connected by trajectories of the geometric flow. Their action is 
simply given by 

Smstanton= ^|AW|, (2.35) 

where AW is the difference of the corresponding values of W at the two critical points. 
Note that AW 7^ 0, in general, since W changes monotonically along the flow lines and 
the instanton action is finite. Then, this yields the standard description of instanton 
solutions of a point particle moving in Euclidean time, but the number of degrees of 
freedom is infinite now, as the evolution takes place in superspace. 

Finally, let us consider the evolution of the volume of spatial slices under the flow. 
In general, it takes the form 

^voI(Md) = \ J d^x^g^idtg.j = / d^x^gg.jE'j (2.36) 

and, therefore, the trace-free part of the driving curvature terms do not contribute to 
the evolution. Otherwise, the volume changes without definite sign. Thus, in princi- 
ple, the volume can bounce along the flow. At the fixed points, where gijE'^^ = 0, the 
volume reaches a local maximum or minimum depending on circumstances. 

As an example, let us consider the Ricci-Cotton flow described by equation ( 2.29[ >. 
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Then, since the Cotton tensor is traceless, the volume of space changes as 



and it can have either sign. Of course, it is possible to normalize this (or any other flow 
that arises in this context) by rescaling the metric with a function of time followed by 
a suitably chosen time reparametrization t — > t(t) so that the volume is preserved 
in time t. This does not resolve the problem, however, since t(t) is not a monotonic 
function of t in general. Thus, the volume does not provide an entropy functional. 

Other entropy functionals might also exist for these flows, generalizing Perelman's 
functional for the Ricci flow [[T5|, but we have not been able to find them. 



3 Bianchi IX model geometry 

All homogeneous space geometries in three dimensions provide consistent ansatz for 
the mini-superspace truncation of the Ricci and Cotton flows and their combination 
thereof. Such spaces follow the Bianchi classification, but for practical reasons we will 
only consider the case of Bianchi IX model geometries that describe homogeneous but 
generally non-isotropic metrics on S^. The corresponding gravitational instantons of 
the four-dimensional Euclidean theory are special in that they admit an SU{2) isome- 
try group and they provide the simplest examples in our study. In this section, we set 
up the notation and present some useful formulas that will enable us to formulate the 
problem as an autonomous system of ordinary differential equations. The Ricci and 
Cotton flows are studied separately here for Bianchi IX geometries. Comparison with 
the gravitational instantons of ordinary gravity will also be made at the appropriate 
places. 



3.1 Some basic facts 

We consider four-dimensional Riemannian manifolds that are foliated by homoge- 
neous three-dimensional spaces of the form 

ds^ = dt^ + Y^j^{t) (cr^y , (3.1) 
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setting N — 1 and N/ = 0. The coefficients jj are taken to depend only on f and cr^ are 
the left-invariant Maurer-Cartan one-forms of 511(2) 

a — sin J? sin ^dcp + cos ^ di? 

£7^ = sin t? cos fdcp — sin i/^ dt? (3.2) 

£7^ = cos I? d^ + dl/' 

with Euler angles ranging as < d < n, < cp < In, < < 4n, which satisfy 

da' + ^e^j^a^ A 0^ = 0. (3.3) 

The three-dimensional leaves are, in general, homogeneous but non-isotropic three- 
spheres. The isometry group is enhanced to 511(2) x ii(l) when any two j/s coincide 
by imposing axial symmetry. Full isotropy requires all 7/s to be equal, in which case 
the symmetry of the model is promoted to 511(2) x 511(2). The volume of the three- 
sphere is 

V = 167tV7i72T3 (3.4) 

and when all coefficients 7, are equal to the volume is expressed in terms of the 
radius LeisV = 2n^L^. 

The Ricci and Cotton tensors are diagonal for all homogeneous geometries and 
this ensures consistency of the reduced models. They take the following form for the 
Bianchi IX class. 



1 

Ru = 



272T3 



Ti 

'11 — 



2(tiT2T3)'/' 



7? - (72 - 73)'J , (3.5) 
7i (271 - 72 - 73) - (72 + 73) (72 - 73)^ / (3.6) 



and similarly for the other two components that follow by cyclic permutation of the 
indices. Also, the Ricci scalar curvature is given by 



R 



27172 + 27273 + 27371-71-72-73] / (3.7) 



2717273 

whereas the trace of the Cotton tensor vanishes, as it can be readily checked. 

With these explanations in mind, we arrive at the following system of ordinary 
differential equations for the metric coefficients 7;(f) 
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as the Bianchi IX mini-superspace model of the combined Ricci-Cotton flow with gen- 
eral couplings. 



3.2 Ricci flow 



When ivqs the Cotton term decouples and one arrives at a variant of the Ricci 

flow as the relevant equation for Hof ava-Lifshitz gravity with anisotropy scaling pa- 
rameter z = 2. In this case, the system becomes 



d7i 
dt 



2A 



2(3A - 1 



A 



w 



3A 



7i 



(3.9) 



and its properties resemble the ordinary Ricci flow on 



d7f 
dt 



-Rnil) ■ 



(3.10) 



Formally, one follows from the other by rescaling the metric with a function of time 
and changing time variable by suitable reparametrization; in such case, the compo- 
nents of the Ricci tensor remain invariant and they assume the same form for the 
rescaled components of the metric. 



The Ricci flow equations (3.101 for homogeneous geometries are well studied in 
the literature following the original work ||37| . For Bianchi IX geometries they take the 
form 



2_d7i 
71 dt 
2_d72 
72 

2 d73 
73 d^ 



1 



717273 
1 

717273 
1 

717273 



(72-73) -7i 



(73-71) -72 



(71-72) -73 



(3.11) 



and coincide with the celebrated Darboux-Halphen system that was introduced by Dar- 
boux in the nineteenth century ||38l and subsequently solved by Halphen i39lj ^ In 
principle, solutions of these equations can be translated into solutions of the original 



system (3.9 1. 



^For the comparison one must consider the variables coi = 7273, C02 = 7173, <^3 = 7i72 and change 
the time coordinate to dT = dt/ 717273- Then, the equations take the equivalent form 



with cyclic permutations for LO2 and 0^3. The identification of the Ricci flow equations for Bianchi IX 
geometry with the Darboux-Halphen system has escaped attention in the mathematics literature. 



18 



It can be shown quite generally that for given initial data 7 ■ the metric will evolve 
towards the configuration 71 = 72 = 73 = by making rounder and rounder until 
the whole space collapses to a point. A particularly simple solution that exhibits this 
behavior is provided by 

li{t) = Hit) = jsit) = \ito - t) (3.12) 

and describes an isotropic metric on whose radius evolves from infinitely large size 
to zero as t varies from — 00 to ig. 

Actually, following the literature ||37ll , the convergence of the flow lines is best 
described in terms of the normalized Ricci flow equation on S'^ 

^ = -Ru + ^Rjr , (3.13) 

which follows directly from equation ( |3.9[ > in the limit A — > — 00, and which can be 
obtained from the ordinary Ricci flow by (yet another) suitable rescaling of the metric 
and time reparametrization. Then, the volume is preserved along the flow and the 
round metric (fully isotropic model with finite radius) arises as fixed point that is 
exponentially reached after infinitely long time, regardless of initial conditions. The 
normalized Ricci flow will also be in focus later in section 4 for different reasons. 

The Darboux-Halphen system is not algebraically integrable when 71 7^ 72 7^ 73- 
All its solutions, however, can be expressed in terms of modular forms (see, for in- 
stance, llQl). The system becomes algebraically integrable when two 7;'s are equal. 
Then, the corresponding three-dimensional space is a three-sphere with axially sym- 



metric metric. Setting 71 = 72 > 73/ the system (3.11) simplifies and exhibits a first 
integral 

111 

(3.14) 



73 7i73 

with arbitrary parameter m. The solution is subsequently described as 

^-^ + arcsinh:^ = = log ^" + - 2m ( ^ ' ) , (3.15) 

m mm 2m — 73 \2m + 73 2m — J 

using another integration constant Iq. The fully isotropic solution with 51/(2) x Sll{2) 
isometry is obtained by taking the limit m ^ 00. 

Remarkably, the same system of equations arises in the description of a class of 
self-dual instantons with SU{2) isometry in ordinary gravity [34J, as well as in the 
description of the moduli space of SU{2) BPS monopoles with magnetic charge 2 [[411 
|42| . In these two cases, which share many features with each other, the ansatz for the 
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Bianchi IX geometries takes the form 



2 2 2 

ds^ = dt^ + a^{t) (^(7^) + h^{t) [a^'^ + c^{t) (cr^) (3.16) 
with a{t) = 7i(0/ ^(0 = 72(0 c{t) = ^3{t) satisfying the same system of equa- 



tions (3.11 1 above. Then, the corresponding axially symmetric gravitational instanton 
with 7i = 72 is the self-dual Taub-NUT metric with NUT parameter m. It can be 
brought into standard form using a radial coordinate r > m with m arcosh (r/m) + 



^/lj^^m){r^^^m) = —t + toso that the solution (3.15 1 becomes 



71 = 72= \/ir-m)(r + m) , 73 = 2m\l'-—^. (3.17) 

r + m 

The fully anisotropic instantons of Einstein gravity with SU{2) isometry correspond 
to the so called Atiyah-Hitchin metric ||4ni42| . which is also the metric on the moduli 
space of charge 2 BPS SU{2) monopoles in general position. These metrics will be 
discussed further in section 6 while comparing instanton solutions of general relativity 
with those of Horava-Lifshitz theory. 



3.3 Cotton flow 



Next, we consider the pure Cotton flow equation^ that arise in the limit — ^ oo- 
They describe solutions of z = 3 Horava-Lifshitz theory when the volume of space is 
very small, i.e., in the deep ultra-violet regime where the Cotton term dominates and 
all subleading relevant operators can be safely dropped from the potential. Then, the 
equations for Bianchi IX geometries take the form 



2wcs d7i 
k'^7i dt 

2wcs d72 
k'^72 dt 

2wcs d73 
k'^73 dt 



1 



(717273)^/' 
1 

(717273)^/' 
1 

(717273)^/' 



7i (271 - 72 - 73) - (72 + 73) (72 - 73)' 



72 (272 - 73 - 7i) - (73 + 7i) (73 - 7h 



73 (273 - 71 - 72) - (71 + 72) (71 - 72)' 



,(3.18) 



and clearly they are much more complicated than the Darboux-Halphen system. It is 
not known whether they are algebraically integrable when 71 7^ 72 7^ 73/ but it can 
be easily shown that, under any initial data 7!*^'', they flow exponentially fast towards 



^The Cotton flow was originally introduced in the literature [36] as an alternative to the Ricci flow 
for studying the existence of constant curvature metrics on three-manifolds. So far it has only been 
applied to homogeneous geometries and its general utility for proving the Poincare conjecture remains 
questionable. 
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the fixed point which is the round metric on and it is conformally flat ||36| . Since the 
Cotton tensor is odd under parity (because of the fully antisymmetric epsilon symbol 
appearing in its definition), there is always an ambiguity in the overall sign of the 
flow equations. Here, we have chosen the sign that takes any metric towards the fixed 
point rather than away from it when Wqs > 0. 

The behavior of the Cotton flow is similar to the normalized Ricci flow, as they 
both preserve the volume of space V = 2n^L?, but the convergence rate is different. 
To compare the two it is sufficient to linearize the corresponding equations around the 
fixed point 

71 = 72 = 73 = Y (^-^^^ 
by considering small perturbations of the metric coefficients 

^i{t) = ^{l + 5x{t)), ^^{t) = ^{l + 5y{t)), (3.20) 
whereas 73 (i) changes accordingly, 

^3{t) = ^{l-5x{t)-5y{t)), (3.21) 

so that the volume of space is preserved. Then, the autonomous system of Cotton flow 
equations ( 3.18[ > becomes to linear order 



d (5x\ _ _ 12k2 n o\ / 
dt\5y)~ ~wcs [o 1 ) \Sy^ 

and the two eigenvalues are equal and negative for cs > 0, 



(3.22) 



1 2k^ 

ensuring stability in all directions around the fixed point. 
The perturbations diminish exponentially fast, as 

Sx{t) = Ae"'/"c , Sy{t) = Be^'^'c , (3.24) 

using arbitrary integration constants A, B and the characteristic time scale of dissipa- 
tion 

Tc = ^ . (3.25) 

Thus, when wqs is very small compared to so that the Cotton tensor dominates 
the flow over the Ricci curvature and the cosmological constant terms, or equivalently 
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when L is very small so that the volume of space is tiny, is very small and the metric 
perturbations dissipate very fast at late time^ 

There is an additional fixed point, which is unique, up to permutation of the axes, 
that arises when 7i = 72 = and 73 = (correlated limit with V held fixed). It corre- 
sponds to a squashed that is completely flattened and has zero Cotton tensor, as it 
can be explicitly checked. Although this configuration is degenerate in one principal 
direction, it has no curvature singularities and it is legitimate to consider. 

Axisymmetric solutions of the Cotton flow can be constructed in closed form, as 
for the Ricci flow. Assuming 

71 = 72 = X— , 73 = ^ , (3.26) 

so that the volume of space is held fixed to V = Itt^I]' , the Cotton flow equations 
( 3.18[ > reduce to a single equation 



which is solved as follows. 



dx 4k^ 1 — 

(3.27) 



-^ + ^* _.3 



— 1 + log 



(3.28) 



with integration constant i^. 

This solution has two branches. For x > 1, x{t) changes from +00 to 1 as i varies 
from —00 to +00; the three-sphere deforms starting from the singular configuration 
7i = 72 = 73 = and gradually becomes rounder until it reaches the isotropic 
fixed point after infinitely long time. For x < 1, x(t) changes from to 1 as ^ varies 
from ^0 = + Tc to +00; in this case, the three-sphere evolves from the singular 
configuration 71 = 72 = 0, 73 = 00 towards the isotropic fixed point. 

Finally, we point out that there is no known solution for the fully anisotropic model 
geometry that is analogous to the general solution of the Darboux-Halphen system. It 
remains an open question whether the system is algebraically integrable and find its 
solution. 



"^We note here for completeness, and it will be used in the next section, that the normalized Ricci 
flow dji/dt = —Rii + Rji/S can be similarly analyzed by considering small perturbations around 
the fully isotropic fixed point. The two eigenvalues also turn out to be equal and negative, but the 
corresponding characteristic time scale for dissipation depends quadratically on L as Tr = /4 with 
respect to the appropriate time coordinate. Comparison with the dissipation rate of the Cotton flow 
will become relevant in section 4. 
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4 Normalized Ricci-Cotton flow 



We will now investigate the combined Ricci-Cotton flow when the effective speed of 
light vanishes by letting A — )■ — oo. In this case, the flow equations take the following 
general form, 

and they become independent of A^y. The driving curvature term is traceless and 
the deformations preserve the volume of space. Thus, the resulting normalized Ricci- 
Cotton flow is a superposition of the Cotton and the normalized Ricci flow with com- 
peting effects, in general, that depend on the sign of zvqs- 



4.1 The general system of Bianchi IX equations 

Using the Bianchi IX ansatz for the three-dimensional geometry, the normalized Ricci- 
Cotton flow equations form an autonomous system of equations for the coefficients 
71/ 72/ 73- Since the volume V = Iti^L^ is conserved, it is convenient to use two 
independent variables x{t) and y{t), 

and also set 

T^^i, F-^. (4.3) 



Then, the general system of Bianchi IX equations takes the form 

+ ( I + ;r) (l + V - 2x^y) ], (4.4) 
+ [^+y)(l + x^y-2xy^)xy\. (4.5) 



It is not known whether this system of equations is rntegrable, in general. How- 
ever, we will be able to characterize its fixed points and study their stability on general 
grounds in order to infer the qualitative behavior of the flow lines. Note that the in- 
dividual components of our system, the normalized Ricci flow and the Cotton flow. 
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can compete with each other. When fi > (wcs > 0)/ they both work in the same 
direction, but when ]i < (wcs < 0) they work against each other and can affect the 
form and stability properties of the fixed points. 

It will also be seen later that these equations can be solved exactly in the axially 
symmetric case x = y. Actually, there are three curves of axial symmetry in the prob- 
lem, but, in practice, it is sufficient to consider only one of these axially symmetric 
cases, since the other two follow by permutation of the principal axes of S^. Thus, 
apart from x = y, we also have x^y = 1 and xy^ = 1, depending on the pair of y/s 
that become equal and reduce the flow equations to a single one. They correspond to 
metrics on with enhanced symmetry SU{2) x U{1) and they all intersect at the fully 
isotropic point x = y = 1. These curves are by themselves flow lines, which, however, 
cannot be crossed by other flow lines; if any two y/s become equal at a given (finite) 
time they will remain equal for ever. Therefore, these three curves provide the barriers 
for six regions in the (x, y) plane where the generic flow lines are confined depending 
on initial conditions. The maximal time range of any given flow line also depends on 
the region in which the flow is confined. Finally, the flow lines along the three curves 
of axial symmetry can reach the fully isotropic point but cannot continue running be- 
yond it. The three curves of axial symmetry are depicted in figure [l| which is restricted 
to the first quadrant of the {x,y) plane so that the metric has physical signature, and 
they intersect at (1, 1). 




Figure 1: The lines of enhanced SU{2) x U{1) symmetry in the (x, y) plane. 
These general qualitative remarks are sufficient to illustrate the evolution of any 
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given initial data on the (x, y) plane, provided that all fixed points are found and the 
arrows of flow lines around them are correctly identified by stability analysis. The 
same remarks apply equally well to the normalized Ricci and Cotton flows which can 
arise separately as special cases. 



4.2 Classification of the fixed points 



We are going now to classify all fixed points of the system ( 4.4[ > and ( 4.5[ > and find the 
critical values of }i that separate their behavior into different phases of stability. 



The isotropic fixed point. The metric on the round sphere is a fixed point of the nor- 
malized Ricci-Cotton flow for all values of ]J. (positive and negative) and corresponds 
to the point 

x^ = y^ = l (4.6) 

on the {x,y) plane. Its stability, however, depends on the values of the parameter }i. 
By considering small perturbations around this fixed point, as 

x{t) = + Sx{t) , y{t) =y^ + 5y{t) , (4.7) 

we find that the linearized system of equations takes the form 



(4.8) 





The two eigenvalues are equal 



= ^2 = - ( 1 + ^ ) (4.9) 



and follow by linear superposition of the corresponding eigenvalues of the normal- 
ized Ricci and Cotton flows discussed in the previous section. 

The fixed point is absolutely stable when }i satisfies the bound 3/}i > —1, i.e., 
^ > or ^ < —3. In these cases, the flow line converge towards the fixed point from 
all directions. Otherwise, for — 3 < < 0, the isotropic point is absolutely unstable 
and the flow lines diverge away from it in all directions. Finally, when }i = —3, the 
eigenvalues are zero and all points in the vicinity of the fixed point are at equilibrium. 
Thus, 



• fi > or fi < —3: absolutely stable fixed point. 
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—3<}i<0: absolutely unstable fixed point. 



These results are in exact agreement with the competition between the Ricci and 
Cotton components of the flow and can be understood by comparing the character- 
istic time scales Tr and Tc that govern metric perturbations at late time^ For fi > 
both components dissipate all metric perturbations exponentially fast, but for }i < 
the Cotton flow contributes differently leading to exponential growth of the pertur- 
bations. For the critical value ^ = —3, the dissipation of the normalized Ricci flow is 
canceled by the exponential growth of the Cotton flow, making zero the characteristic 
matrix of metric perturbations. 



Anisotropic fixed points. For negative }i there are additional fixed points that corre- 
spond to particular axisymmetric metrics. As such, they appear in three copies related 
by permutation of the axes of and they are located on the lines of axial symmetry 
in the (x, y) plane. In particular, they arise 



• on the diagonal x = y with = a/— 3/ }i , 

• on the branch x^y = 1 with x^ = a/— 3/ }i , 

• on the branch xy^ = 1 with y^ = a/— 3/ }i . 

Note that all these points coalesce with the fully isotropic point when ji assumes the 
critical value }i = —3. 

It suffices to perform stability analysis around one of these fixed points, say the 
one located on the diagonal line, since the results will be identical for all of them by 
the symmetry of the problem. Using small fluctuations around the fixed point 




x* = y*=y--, (4.10) 
so that x{t) = X* + 5x{t) and y{t) = y* + 5y{t), the linearized system takes the form. 




^According to definitions we have the following relations 

^ = ^='tk ^^^^ = '^ = ^- 

The characteristic time scales tr and refer to the original time coordinate t but their ratio is the same 
in the time coordinate t. 
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where 



^1 

C2 





(-1 


3 




fi 


'(-- 


3 





3/2 



3/2 



3/2 



(4.12) 
(4.13) 



The two eigenvalues are and ^2 and they are unequal offering various possibili- 
ties in general. They both vanish for fi = —3. Otherwise, we have the following cases 
depending on the sign of the eigenvalues: 

• —3<}i<0: saddle fixed point with ^1 < < ^2 / 

• —6 \/2 < }i < —3: saddle fixed point with ^2 < < ^1 , 

• }i < —6 yjT:. absolutely unstable fixed point with > . 

Thus, the axisymmetric fixed points are never absolutely stable. They are saddle or 
unstable depending on ]i. 



Totally anisotropic fixed points. It is not obvious from the beginning whether there 
are any fixed points with 71 7^ 72 7^ 73- Close inspection of the equations, assisted by 
numerical scanning, reveals the presence of two totally anisotropic fixed points that 
coexist with the axially symmetric anisotropic fixed poinlj^when ji < —6 \/2. In fact, 
by the Z3 symmetry of the problem, there are six such additional fixed points, but 
we only focus attention on two of them appearing symmetrically left and right of the 
diagonal line x = y in the lower two (out of the six disconnected) regions shown in 
figure 1; their presence should not be confused with the mirror images of the axially 
symmetric anisotropic fixed point discussed earlier. 

The characteristic property of the totally anisotropic fixed points is that their Ricci 
scalar curvature vanishes and their location on the (x, y) plane is given by 

x+-i= = -|=y+-L with x^y. (4.14) 

Of course, one can always find Bianchi IX metrics with zero scalar curvature by im- 
posing the appropriate algebraic condition on the metric coefficients, but these are not 
fixed points of the flow lines for general values of ji. Remarkably, they are real solu- 
tions of dx/dr = = dy/dr with x ^ y, which coexist with the axially symmetric 

^We thank Christos Sourdis for pointing out the presence of these additional fixed points and thor- 
oughly investigating their properties. This analysis was missed in a previous version of our paper and 
we are indebted to him for providing all the details. Similar results also appeared in li43J . 
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anisotropic fixed point when ji < —6 ^/l. They are not present when ft > —6 
since there are no real solutions to equation ( |4.14[ > in that case. We also note that when 
}i = —6 \/2 the totally anisotropic fixed points coalesce with the axially symmetric 
anisotropic fixed point. Furthermore, there are no other fixed points in the problem. 

The location of these fixed points is depicted in figure |2j Here, we plot the ratio 
x/y of the anisotropic fixed points as function of }i (actually The horizontal line 
represents the axially symmetric anisotropic fixed point that exists below and above 
the critical value of }i. 
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Figure 2: Relative location of all anisotropic fixed points for ^ < 0. 

At this point, it is instructive to consider the Ricci scalar curvature of the axially 
symmetric anisotropic fixed point = = ^—3/ }i, which turns out to be 

« = -|?('' + ^^\/-f)- ("^> 

When ji = this point is at infinity and the curvature is zero. As ji varies from 
to —3 the axially symmetric anisotropic fixed point is approaching the isotropic fixed 
point and the curvature increases monotonically. The curvature becomes maximal 
when these two points coincide at the critical value }i = —3, and, then it decreases 
monotonically as }i varies from —3 to —6 \/2. It becomes zero at the other critical 
value ji = —6 and, then, it turns negative for fi < —6 \/2. The value ji = —6 yjT. 
is also critical for the creation of the totally anisotropic fixed points, which pop out 
symmetrically from the diagonal and have zero curvature for all lower values of the 
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parameter fi. 

These additional fixed points appear to be saddle points, as can be verified by nu- 
merical investigation for different values of ^. It is not easy to obtain closed formulas 
for the eigenvalues of the characteristic matrix describing small perturbations around 
them. However, their stability properties are important for constructing instanton so- 
lutions of Horava-Lifshitz gravity, as will be seen in detail later in section 6. Another 
characteristic property of these points that will also be used later is the universal value 
of the gravitational Chem-Simons functional. Explicit calculation shows that W^s for 
the fully anisotropic fixed points is 

Wcs = — , (4.16) 

which is independent of ^! It actually coincides with the value of W^s for the axially 
symmetric anisotropic fixed point when ji = —6 yjT. and it is ten times larger than the 
value of Wcs evaluated at the totally isotropic fixed point. 

An important remark is in order at this point. The equations that determine the 
fixed points of the normalized Ricci-Cotton flow provide only the traceless part of the 
classical equations of motion of topologically massive gravity, leaving the trace unde- 
termined. Then, depending on the value of their Ricci scalar curvature, these points 
also satisfy the trace equation R = 6Aw for appropriately chosen effective value of 
Aw- Fixed points with positive, negative or zero Ricci scalar curvature are vacua of 
topologically massive gravity with positive, negative or zero cosmological constant, 
respectively. As a result, the fixed points of the Ricci-Cotton flow with general cou- 
plings, which are the vacua of topologically massive gravity, are expected to be less 
than the fixed points of the normalized flow, and, in fact, they can be obtained from 
them in certain ways This, however, does not make our analysis in section 5 redun- 
dant since their location, volume and stability properties also depend crucially on the 
flow equations we are considering in each case. 

Summarizing the results for the normalized Ricci-Cotton flow, we note that the 
fully isotropic metric on is the unique fixed point for ji > that attracts all trajec- 
tories starting from any point in the first quadrant of the {x,y) plane. For ji < 0, there 
are various possibilities that result to attractive or repelling directions around the fixed 
points. Note that the isotropic point becomes absolutely unstable when —3<}i<0, 
in which case the anisotropic fixed point is a saddle that attracts partially the flow 
lines. For —6 \/2 < }i < there are four in total fixed points, including their Z3 
mirrors, whereas for fi < —6\/2 the total number of fixed points is ten. 



^"^The mathematical problem one has to solve to determine the set of fixed points of massive gravity 
for given Aw is to fix R instead of the volume of the normalized fixed points, as functions of }i and L, 
and deduce from it the allowed range of }i for the selected set of fewer points. 
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4.3 Phase portraits of the flow 



The qualitative behavior of the flow lines is illustrated by the phase portraits shown 
below for all possible values of fi. 

For positive }i, which is qualitatively the same as for the normalized Ricci flow and 
the Cotton flow, the phase portrait is given in figure |3] 
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Figure 3: The flow lines for ^ > 0; here, }i = 2. 

For negative }i the stability properties of the fixed points is different in the intervals 
—3 < }i < 0, —6\/2 < }i < —3 and }i < —6 \f7. and all these possibilities are repre- 
sented in figures |4| |5] and [6| respectively. For later reference, it is important to realize 
the existence of trajectories interpolating between different fixed points. In figure |4] 
there is only one such flow line since one fixed point is unstable and the other a sad- 
dle. The same applies to figure |5] where the interpolating flow line connects a stable 
fixed point with a saddle. The picture changes drastically in figure [6] since there are 
infinitely many flow lines interpolating between an unstable and a stable fixed point. 
Also, in this case, there are additional fixed points away from the diagonal, which are 
saddle points. Figure [6] contains two such totally anisotropic fixed points located at 
X ~ 0.19, y ~ 1.75 (and x ~ 1.75, y ~ 0.19) for }i = -10 < -6^/2c^ -7.56. There are 
flow lines connecting these fixed points with the other two lying on the diagonal, but 
they are not easily seen on the phase portrait due to numerical deficiency. 
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Figure 4: The flow lines for — 3 < ^ < 0; here, ^ = —2. 
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Figure 5: The flow lines for —^\j7.<}i< —3; here, }i 
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Figure 6: The flow lines for }i < —6 -v^2; here, }i = —10. 
4.4 Axisymmetric solutions 

Here, we present the exact axially symmetric solution of the normalized Ricci-Cotton 
flow, setting x = y. As such, it generalizes the axially symmetric solutions of the 
normalized Ricci and the Cotton flows discussed in the previous section. The system 



(4.4 1 and (4.5 1 reduces to a single equation 



^"-'^^ + (4.17) 



dr fi \x^ / V^^ 3 

and the flow takes place on the diagonal line of the {x,y) plane. The flow connects 
either the origin or infinity with a fixed point or it can extend between two different 
fixed points. These regions do not communicate with each other and the time interval 
that supports the solutions depends on ji and the choice of trajectory. 

The behavior around x = and x — ^ oo is universal and can be extracted directly 



from equation (4.171, 



X 



/ , N 1/6 

(^^(T-To)j , asx^O, 



X — 5T, as X — > 00. 
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Around the fixed points the time dependence is exponential and determined by the 
eigenvalues ^ as x — ^ exp(^T). The relevant eigenvalue for the isotropic point is 
given by (4.9) and for the anisotropic by ( |4.12[ ). 



Case I: ^ > 0. The solution behaves similarly for all positive values of }i, but it looks 
different on the two sides of the isotropic fixed point = 1. We find that 



X > 1 



-00 < T < +00 as X decreases from +oo to 1 , 



X < 1 : Tq < T < +00 as x increases from x(to) = to 1 



The exact solution reads 



-3x + „ . ^ log f X^ + X + 1 ) + 2 



2 (^2 -3^ + 9) 



81 



log I X — 1 1 + 
3 



2(^3 + 27) 



log (J^x^ + 3^ 



arctan 



2x + l 



, . I arctan x , , 
11^ + 27^11 V V 3 



(4.18) 



and it is represented by figure |7] with appropriately chosen integration constant t^. 




Figure 7: The dependence t(x) for positive ji; here, ji = 1. 
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Case II: }i < 0. In this case the solution depends on the particular value of }i. First, 
we present the result for the critical value }i = —3, which is simpler. 



T- 



-3x + 



2(x 



1 2x + l 7, , 

— + — arctan — — log x - 1 

1) a/3 a/3 4 ^' 



3 1 
+ -log(x + l) + -log + x + 1 



(4.19) 



Figure 8: The dependence t(x) for the critical value critical ji 



All fixed points coalesce to the isotropic and the eigenvalues vanish, so that there 
are no arrows infinitesimally close to this point. However, this behavior is lifted at 
second order in perturbation theory and there are arrows pointing from large to small 
values of t. This particular case is depicted in figure |8| which also shows the range of 
T in the two branches. 



For ^7^—3 the expression becomes much more involved and reads 

i_ , ^ , iA , V3}1^ 



-3x + 

}i + 3 
81 



l{}fi-3}i + 9) 
log I X — 1 1 + 



log X^ + X + 1 + 



arctan 



2 x + l 



27fi 



2(^3 + 27) 



log fix + 3 



2(^3 + 27) 



3 



log 



X 



log 



1 + X 



(4.20) 



Then, depending on whether — 3< jd < Oor ji < —3, the function t(x) looks different. 
In all these cases the solution has three branches but the range of time is not the same. 
We have, in particular, the following behavior depending on ji: 
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For —3<}i<0 the two fixed points are ordered as x^^° = 1 < xf^^° = a/— 
with the isotropic being repulsive and the anisotropic attractive. The exact solu- 
tion is represented by figure |9| 




Figure 9: The dependence t(x) for — 3 < ^ < 0; here, y. = —1. 



For < —3 the two fixed points are ordered differently, as x 



aniso 



-3/}i < 



1. The isotropic point is now attractive whereas the anisotropic is repul- 



sive. Then, the exact solution along the diagonal is represented by figure 10 and 
there is no essential difference for }i below or above the value —6 \/2 . 




Figure 10: The dependence t(x) for }i < —3; here, ji = —5. 

Note that the three branches shown in figures 9 and 10 range differently: the first 
has support on a semi-infinite time interval, whereas the second and third to the 
right are eternal solutions that exist for all time — oo < t < +oo. (Despite ap- 
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pearances, caused by numerical deficiency, the spike at x = 1 extends to infinity 
in figure 9). 



5 Ricci-Cotton flow with general couplings 

Let us now examine the Ricci-Cotton flow equations for Bianchi IX model geometries 
with arbitrary couplings by letting A and A^v take arbitrary values. The system of 
equations that needs to be studied is provided by ( |3.8[ > 

with 

T=4^ F='^- (5-2) 

The definition of T and }i resembles that for the normalized flow in section 4, but it 
does not include the rescaling by the characteristic length of space. 

Since A < 1/3, we will confine our discussion to the case of non-negative cosmo- 
logical constant. Aw > 0, so that the effective speed of light is real, and investigate the 
structure of the fixed points and their stability properties. The choices A^ > and 
Aw = will be discussed separately, although the latter can be obtained as limiting 
case of the former. Mathematically it is also interesting to consider the more general 
situation, without imposing any restrictions on A and Aw, but these cases will not be 
included herj^ 



5.1 The axisymmetric Bianchi IX model 

The flow of the metric coefficients yi{t) does not preserve the volume of S"^ in this case, 
and, therefore, the three coupled equations are rather difficult to investigate in all gen- 
erality with 7i 7^ 72 7^ 73- Restricting attention to axially symmetric configurations 
simplifies matters without shadowing too much the rich structure of the system. Our 
analysis will be based on the qualitative theory of dynamical systems, as in previous 
sections, but because of the mathematical complexity of the equations it is not possi- 
ble to obtain explicit solutions in closed form, apart from the isotropic solution. Also, 
the reader should be aware of the limitations: the fixed points can become unstable 
in other directions, when axial symmetry is relaxed, and the conclusions drawn here 

^^The nature of the corresponding fixed points changes drastically compared to the cases that are 
discussed in this paper. However, when A^y < it is more appropriate to consider Bianchi type VIII 
model geometries rather than Bianchi IX. 
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may be altered and not be as general. Of course, this is part of a more general crit- 
icism for mini-superspace models compared to metric deformations with an infinite 
number of moduli taking place in the entire superspace. 

With these explanations in mind, we are going to study the system of two equa- 
tions 



= _^(,,_,,)_Jii±^2| + ?i:ill + ^, ,5.3) 

71 dr m2 4(3A — 1) 7| 3a — 1 72 3A — 1 

72 dr 2^ 73 ^ ^ 4(3A - 1) 72 3A - 1 72 3A - 1 ' ^ ' 

which is obtained from ( |5.1[ > by setting 72 = 73. The existence and properties of the 
fixed points depends crucially on the values of as in other examples considered so 
far. In particular, for ^ < 0, an anisotropic fixed point will coexist with the isotropic 
one. 

There is an exact solution of these equations which is available for all Aw > and 
describes the evolution of the isotropic metric on S^. This possibility does not arise 
for the normalized Ricci-Cotton flow, since the isotropic metric is a fixed point. In 
particular, setting all 7, = 7, we have the following result: 

7W=Aexp(3^x)+ji- forA„>0 (5.5) 

and 

7(0 = -^(^;^I:^(t-to) for Aw = (5.6) 

with A and Tq being arbitrary integration constant. Since A < 1/3, the metric flows to 
the isotropic fixed point (to be discussed next in detail for more general trajectories) 
as T — > +00. 

The difference between Aw > and Aw = is reflected in the life-time of the 
solutions. When Aw > the solution has two branches: on the first branch A > 
and the solution is eternal existing for all —00 < t < +00; it describes a round sphere 
with infinite radius in the infinite past flowing towards a round sphere with radius 
set by Aw in the infinite future. The second branch corresponds to A < in which 
case the solution exists for < t < +oq, with appropriately chosen t^, so that the 
sphere starts from zero radius and reaches the fixed point as f — > +00. On the other 
hand, when Aw = 0, there is only one branch as the solution exists for tq < t < +00, 
interpolating between a fully collapsed configuration at t = tq to a sphere of infinite 
radius in the infinite future. 
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5.2 Classification of the fixed points 



First, we consider the case of non-vanishing cosmological constant and reserve the last 
subsection to study Aw = separately. 



The isotropic fixed point. For > 0, there is a natural length scale in the problem 
that gives rise to the isotropic fixed point of the flow, irrespective of the sign of }i, 



4A 



(5.7) 



This follows easily from the system of equations without assuming any restrictions on 
7i and 72; it can also be shown that it is a fixed point of the more general system of 
equations with 7i 7^ 72 7^ 73- 

Linearizing around it as 



71 (0 



(l+Sxit)), 72(f) 



(i + ^y(0). 



4Aiv 4Aiv 
we find that the small perturbations satisfy the characteristic matrix equation 



(5.8) 



d_ 
dr 



( 8 r— 8A-3 8 r— ^4A-1 \ 



= A 



The eigenvalues are 



4: r— 4A-1 



3A-1 

4 r— 2A-1 



/Sx\ 



(5.9) 



A 



w 



3A-1' 



^2 = -4A 



(1 + ^^) 



(5.10) 



and so ^1 < whereas ^2 can take all values, positive or negative, depending on }i. 

Keeping Aw fixed and varying ji we obtain the following characterization of the 
isotropic fixed point: 



^ > or ^ < — 3\/Aw: absolutely stable, 
— 3v^Aw < }i <0: saddle point. 



Note the emergence of a critical value, ^ = — 3i/Aw, where the two eigenvalues ^1 
and ^2 vanish, separating stability from instability along the corresponding eigen- 
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directions Also note for completeness that if we were allowing A > 1/3, the 
isotropic point would never be absolutely stable (it would be absolutely unstable or a 
saddle point in the respective intervals of the ^-line.). 

Actually, one can go further and investigate whether the exact isotropic running 
solution j{t), given by (5^, is stable against small fluctuations. 



7/(0 =7(0+^7.(0 



(5.11) 



acting as attractor of nearby trajectories. Thus, given a small tube around the trajectory 
7(i), one is interested to know if any other trajectory with initial conditions inside this 
tube will remain there after some time and what is the size of tube that guarantees this 
attractor property. The differential equations for 5^i{t) are most conveniently stated 



using 7(f) rather that t as flow time. Then, within the axially symmetric ansatz (5.3 1, 
the linearized equations take the following form. 



d 

d^ 



-^7i(0 
-^72(0 



2((57i - 5^2) 
7(4Aw7 - 1' 

h\ - hi 

7(4Aw7 - 1) 



2(3A 



H\/7 
2(3A-1) 
^^7 



+ 4A-1 



+ 



+ 4A- 1 



+ 



4Aw^7i 
4Aw7 

4Aw^72 
4Aw7 - 1 



1 



(5.12) 



(5.13) 



generalizing the characteristic matrix equations of small perturbations around the 
isotropic fixed point. These equations apply for all Aw including the special case 
Aw = that will be discussed separately. 

Solutions of this system can be expressed as functions of t through 7(f). Since ^7/ 
are required to be small, for validity of the linearized analysis, the attractor property of 
the isotropic trajectory appears to be very limited. This behavior can be seen schemat- 
ically in the phase portraits of the flow that will appear in the next two subsections. 



Anisotropic fixed point. When }i < 0, there is an additional fixed point associated 
with the axially symmetric metric with coefficients 



36^2 

{pi^ + 27Awf ^^ + 27Aw" 



71 = — — — ^2' 72 = 73 = ..7 , • (5-14) 



There are no other restrictions on the values of ^ for the existence of this second fixed 
point. 

^^The critical value of /i is similar to that found for the normalized Ricci-Cotton flow; direct com- 
parison can be made by replacing ^/A-y^ with 1 / L and rescaling ja with the characteristic length of 
space. In both cases, the critical value of f/ occurs when the competing effects of the Ricci and Cotton 
deformations are balanced exactly. 
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Notice that this new fixed point and the isotropic one will coalesce if ^ = — Sa/Avv. 
For — 3\/Aw < }i < 0, the anisotropic fixed point has 71 < 72, whereas for }i < 
— 3\/Aw it has 71 > 72. Thus, — 3\/Aw appears as a critical value of 

By considering small perturbations around the anisotropic fixed point, as 



7i(0 



+ 27 Aw] 



{1 + Sx{t)), 72(0 



^2 + 27A 



w 



:i + Sy{t)), (5.15) 



we find the characteristic matrix of the linearized system with respect to t 

(9A - 2)}i^ - 27Aw(3A - 1) -(18A - 5)^^ ^ 27Aw(6A - i; 



27(3A-1) I -i(9A-5)^2^^Aw(3A-l) (9A - 4)/^2 _ 27a^(3A - 2) 



The corresponding eigenvalues are 

1 



18(3A- 1^ 



2(3A - 1)^2 _ 27Aw(2A - 1) ± Va 



where 



A = 6(3A-1)(2A-1)^^-72AwH^(3A- 1)2 + 243A5v (i2A.^-6\ + 1 



(5.16) 



(5.17) 



(5.18) 



The eigenvalues ^± are real^ Since their product is given by 



(f/2 - 9Aw) (f/2 + 27Aw) 
162(3A-1) 



(5.19) 



we note the appearance of a critical value jd, which is the same as for the isotropic fixed 
point, }i = —3y/Aw. Then, for ji < —3y/Aw the anisotropic point is saddle. On the 
other hand, in order to examine the stability of this fixed point for — 3\/Aw < }i < 0, 
we consider the sum of the two eigenvalues. 



9(3A-i; 



2(3A-l)//2_27A(2A-i; 



(5.20) 



which is now negative. Therefore, for — 3v^Aw < ^ < 0, the anisotropic fixed point is 
absolutely stable. Summarizing all results obtained above, we have the following: 

^^It follows by noting that the two roots of A occur at 

3Aw 



A 



2(3A-1)(2A-1) 



4(3A-1)^± ./2(3A-1) 



and they are complex for A < 1/3. Thus, A has the same sign as the coefficient of its /;^-term, which is 
positive. 
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• An isotropic fixed point exists for all fi and it is absolute stable when ^ > or 
}i < — 3\/Avv- For —3y/ Aw < < it is a saddle point. 

• An anisotropic fixed point exists for all ^ < 0. It is absolutely stable for —3y/ Aw < 
}i < and saddle for ji < —3^/Aw, which is reverse to the behavior of the 
isotropic fixed point. 

As can be seen there are similarities as well as some differences with the classifica- 
tion of fixed points of the normalized Ricci-Cotton flow. 



5.3 Phase portraits of the flow 

A qualitative picture of the flow lines is provided by three consecutive phase portraits 
for different values of the parameter ji. In all drawings we choose Aw = 0.25 and 
A = 0.1, and so the three regimes are ^ > or —1.5 < ^ < or ^ < —1.5. The 
isotropic fixed point appears at 7i = 72 = l/4Aw = 1- 

First, we consider the case ji > that exhibits only one (isotropic) fixed point, as 



shown in figure 11 



1.5 



0.5 











r r > 



1 r ' 



' 1 r 



* •> 1 
V \ a' 



0.5 



1 

71 



1.5 



Figure 11: The flow lines for fi > 0; here, fi = 1. 
Next, we consider the case < that exhibits an additional (anisotropic) fixed 



point and make the following choices for the plots shown in figures 12 and 13 
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Figure 12: The flow lines for — 1.5 = — 3\/Avv < ^ < 0; here, — 
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Figure 13: The flow lines for ]i < —3\/Av^ = —1.5; here, ji = —2. 
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fi = —1: the anisotropic fixed point occurs at 71 ~ 0.60, 72 = 73 ~ 1.16. 
}i = —2: the anisotropic fixed point occurs at 71 ~ 1.24, 72 = 73 ~ 0.84. 



5.4 The special case Aw = 

Setting Aw = corresponds to taking the effective speed of light equal to zero, while 
keeping A < 1 /3 arbitrary. It should be contrasted with the normalized Ricci-Cotton 
flow, which is independent of Apy and has also zero effective speed of light. 

The isotropic fixed point is now pushed to infinity and corresponds to a round 
with infinite radius. This is also apparent from the exact isotropic solution ( |5.6[ > that 
converges to it after infinitely long time. However, it is not strictly speaking a fixed 
point of the flow lines, since d7i / d^ and d72/df do not vanish there when Aw = 0. 

When < there is an anisotropic fixed point of the axially symmetric flow for 



36 , 

71 = ^ = 472 = 473 . 



(5.21) 



This is always a saddle point because the corresponding eigenvalues of the character- 
istic matrix of small perturbations are real for A < 1 /3 and their product is negative. 



The results follow setting Aw = in the expressions we had before (see e.g. (5.19 1). 



2 
1.5 
1 

0.5 








4 i 



/ i k V V 



0.5 



1 



1.5 



Figure 14: The flow lines for pi > 0; here, fi = 1. 
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Figure 15: The flow lines for fi < 0; here, }i = —5. 

We include two phase portraits of the flow lines that are characteristic for ^ > 
and ji < 0, respectively, choosing A = 0.1. The case ^ > is shown first in figure 



14 Next, we consider the case }i < that exhibits an additional (anisotropic) fixed 



point and make the choice fi = —5. The anisotropic fixed point occurs at 71 ~ 1.44, 



72 = 73 — 0.36 as shown in figure 15 



Finally, concluding this section, we end up with an interesting observation that 
arose in our study of the flow lines. Along these lines, the volume of space changes, 
but the dependence on t is not monotonic in general. It is therefore interesting to 
inquire in this context for the existence of bouncing solutions for which the volume 
reaches a minimum and then increases in time. Although this behavior is not generic, 
it seems to arise along particular flow lines that can be found by numerical scanning. 



An example of this kind is provided in figure 16 for an appropriate choice of initial 
data and couplings. Cases with full anisotropy, but with initial conditions close to 
axial symmetry, also seem to lead to bounces, including minima with very small vol- 
ume. Similar conclusions hold when Aw > 0, with or without axial symmetry, but we 
have not been able to obtain any quantitative characterization of the phenomenon so 
far. 

The bouncing solutions can be regarded (in some sense) as the Euclidean space 
analogue of bouncing models in standard cosmology, which provide a viable alterna- 



44 



> 




-5 -4.5 -4 -3.5 -3 -2.5 -2 -1 .5 -1 -0.5 
Log(T) 

Figure 16: Bouncing solution with initial data = 28, 72*^'* = 73*^^ = 4 for A = 0.26, 
Aw = and }i = —5. 



tive to inflation. Matter bounces have already appeared in studies of Horava-Lifshitz 
cosmology (see, for instance, ||44ll45ll46| ), but they are also non-generic. It remains to 
be seen whether they have any special meaning and implications for the models we 
study here, although their occurrence does not require any matter couplings, as they 
are purely geometric, and, hence, different from those arising in cosmology. 



6 Space-time interpretation of the flow lines 

The solutions of geometric flows are shown as flow lines in the various phase por- 
traits that have been drawn. Explicit solutions were obtained in special cases, whereas 
more general solutions can only be described pictorially. The problems that will be 
addressed in this section are the selection of flow lines that can qualify as regular 
gravitational-instanton solutions in four dimensions and the completeness of the cor- 
responding space-time metrics. It will also be useful in this context to compare SU{2) 
instanton solutions of Horava-Lifshitz theory with those of ordinary Einstein gravity. 

According to the analysis of section 2.3[ only the flow lines that interpolate be- 



tween two fixed points (when more than one fixed point is present in our models) 
qualify as instantons. They are indeed finite-action solutions, and this property is 
sufficient to determine the global structure of space-time and its asymptotic behav- 
ior, and render the corresponding metrics complete. Our main result in a nut-shell is 
that all gravitational instantons with SU{2) isometry are globally R x describing 
a smooth deformation of as t runs from —00 to +00, without ever encountering a 
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singularity. The details are given below together with the complete classification of grav- 
itational instantons with SU{2) isometry in Horava-Lifshitz gravity with anisotropy 
parameter z = 3. We will also compute their Euclidean action and determine their 
moduli spaces. In most cases we have explicit solutions. There are also a few other 
solutions that are shown to exist, but we have not (yet) been able to obtain expressions 
for their metric in closed form. 

The variant of the theory with anisotropy parameter z = 2 will not be especially 
discussed, since it is clear that it exhibits no instanton solutions (of the type we are 
considering here) with SU{2) isometry. Recall in this case that the relevant equations 



are provided by the Ricci flow on S^, which takes the form (3.9 1 in proper time with 
arbitrary parameters A < 1 /3 and A^y > 0. Its fixed points are determined by the 
equations 

and clearly there is a unique solution given by the constant-curvature metric with 
R = 6Aw- The absence of other fixed points, which is also implied on more gen- 
eral grounds by Poincare's conjecture for S^, shows that there can be no finite-action 
instanton solutions in this case. 

Thus, in the following, we focus on instanton solutions of the z = 3 theory and 
explain their properties, as outlined above. 



6.1 Global structure and completeness of the metrics 

First of all we examine the occurrence of singularities that can appear at finite proper 
time t (or t, since the two are simply related by rescaling) and render the Euclidean 
space-time manifold incomplete. Singularities arise when some or all of the metric co- 
efficients of the Bianchi IX model geometry vanish and they are classified, in general, 
as nuts and bolts. Such singularities are intimately related to the fixed points of Killing 
vector fields by geometrical reasoning and they are independent of the gravitational 
equations. Following ||47l |33ll we recall that the structure of the fixed point set of a 
Killing vector field acting on any four-dimensional Riemannian manifold with met- 
ric is determined by the rank of the 4x4 matrix Vpi^y. This is an anti-symmetric 
matrix (since its symmetric part vanishes identically by definition of a Killing vec- 
tor) which can have rank 4 or 2; rank is excluded for, otherwise, the vector field 
vanishes everywhere. In the former case there are no directions left invariant at the 
tangent space of the fixed point, which, thus, appears to be isolated and it is called 
nut. In the latter case only a two-dimensional subspace of the tangent space at the 
fixed point remains invariant under the action of the Killing vector field, whereas the 
two-dimensional orthogonal complement rotates into itself. Then, the fixed point set 
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is provided by this invariant two-dimensional subspace and it is naturally called bolt 
(it is typically a two-sphere, as in Bianchi IX geometries). 

Nuts and bolts lead to incomplete manifolds, in general, but in certain cases the ap- 
parent singularities can be removed and provide regular and complete metrics with no 
curvature singularities. It all depends on the form of the metric as these singularities 
are approached. A removable nut singularity contributes one unit to the Euler num- 
ber X of the four-manifold and a removable bolt singularity contributes two units 
||47|. following the theorems on fixed points. This counting applies to compact four- 
manifolds without boundary, but it also generalizes to non-compact spaces when the 
Killing vector field is either everywhere tangential to the boundary (as in space-times 
with homogeneous spatial sections that we are considering here) or is everywhere 
transverse. Thus, if 7^ 0, any Killing vector field will have at least one fixed point. 
No fixed points imply that = 0. Instanton solutions of Einstein and Horava-Lifshitz 
gravity are quite different in this respect, as will be seen shortly, having important 
implications for their global topological structure. 

Let us briefly review when such singularities can be removed from a Riemannian 
four-manifold without referring to any specific theory or any solutions at the moment. 
Using locally the Bianchi IX form of the metric ( 3.16[ > 



2 2 2 

ds2 = di^ + a^{t) (cri) + + c^[t) , (6.2) 

we suppose that a singularity (nut or bolt) occurs at some finite proper distance, say 
i = 0. It is well known that the metric has a removable nut singularity provided that 
near i = all metric coefficients vanish as 

a^{t) = b^{t) = c^{t) = as t^O. (6.3) 

In this case we have a coordinate singularity of the polar coordinate system in 
centered att = 0, which is simply removed by changing to a local Cartesian coordinate 
system near the point t = and adding it to the manifold. Also, it is well known that 
the metric has a removable holt singularity provided that near i = two of the metric 
coefficients (say a^ and iP-) become equal and the third vanishes as 

a^[t) = b^{t) = finite , c^{t) = -n^t^ as t ^0 with n G Z . (6.4) 

Then, a^ = b^ = implies a^{o-'^)^ + b^{cr^)^ = R§(di?2 + sin^Mtp^), which is the 
canonical metric, while the dt^ + c^{a^)^ part of the four-dimensional metric be- 
comes dt^ + {n^t^ / 4:)dxp^ near t = 0, keeping and cp constant. In this case, the 
topology of the manifold is locally IR^ x and the factor shrinks to a point on 



4:7 



as i — > 0. By adjusting the range of so that nip/ 2 runs from Q to In, the apparent 
singularity at f = becomes a coordinate singularity of the polar system in and 
can be removed as before. In all other cases the singularities cannot be removed and 
the manifold is incomplete. 

The above reasoning is purely geometrical without reference to any field equations. 
Thus, different gravitational theories for Euclidean space-times of the form ( |6.2[ > may 
or may not lead to removable singularities at the fixed points of a Killing vector field. 
This depends on the way that the metric coefficients approach zero in the vicinity of 
a singularity and it is sensitive to the dynamics. The space-time singularities of Eu- 
clidean Einstein and Hof ava-Lif shitz gravity (if they are present) will follow different 
power-law behavior, which, in turn, will affect the completeness of the corresponding 
metrics. Thus, the absence of non-removable singularities provides a natural selection 
for the physically admissible solutions in those theories. 

In Euclidean Hof ava-Lifshitz gravity, a singularity can only arise if an eligible flow 
line reaches the boundaries - including the origin - of the physical parameter space, 
namely the two wedges of the first quadrant in the {x,y) (or (71,72)) plane used in 
the phase portraits. Then, the flow becomes extinct as it cannot be continued beyond 
that point. Such singularities, if they are present, will arise at finite proper time (say 
t = tQ, but we can always set tQ = without loss of generality). This is obviously 
so because such singular points can also act as initial data for the time-reversed flow 



at a given finite instance of (proper) timd^ Such possibilities should be ruled out 



by the theory, unless the singularities are removable nuts or bolts, for, otherwise, the 



space-time metric will be singular. Using the Bianchi IX form of the metric (6.2 1 with 



a 



^ — 7i, = 72 and = 73 (like in (3.1)), which is the appropriate choice in this 



case, we may set (as in nuts) 

7i = /3ifPS 72 = /32^PS l3 = ht^' (6.5) 
or alternatively (as in bolts) 

7i = i6i, 72 = 162, l3 = ht^' (6.6) 

and determine the allowed values of the coefficients /3, and the exponents p, as f — 0. 
In all cases we find that the first-order system of Bianchi IX equations for Hofava- 
Lifshitz gravity with general couplings does not lead to removable nuts or bolts. Only 
non-removable singularities can arise along the flow lines and they seem trouble- 



^"^Proving finiteness of extinction time for the solutions of geometric flow on certain three-manifolds 
is an intricate mathematical problem that will not be addressed in all generality, since we are only con- 
sidering homogeneous geometries on S^. It is a key point in Perelman's proof of the Poincare conjecture 
based on Ricci flow 1151 and it is not yet clear how it may generalize to the Ricci-Cotton flow. 
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Recall, however, that instanton solutions are rather restrictive, since they are 
described only by those trajectories that interpolate between two fixed points. There- 
fore, their metrics would be singular if any one of the fixed points were singular. But 
this is a contradiction of terms and cannot happen, since a fixed point, unlike a singu- 
larity, is only reached at infinite proper time, f — > ±00 (with sign that depends on the 
direction of the flow), and not at at finite time Thus, the instanton solutions protect 
themselves from the singularities that may otherwise arise by moving along generic 
flow lines. Whenever instanton solutions exist, their spaces will be always complete 
without any singularities. 

In Euclidean Einstein gravity nuts and bolts are important elements in the theory 
of gravitational instantons, since all known solutions exist thanks to their presence. 
In this case, there is a certain class of gravitational instantons ( 6.2[ > that follows from 
the Ricci flow equations ( 3.11[ > in proper time t with 71 = fl, 72 = & and 73 = c, as 
explained in section 3.2. They include the trivial flat-space metric associated with the 
isotropic solution (3.121, having fl^ = fc^ = = everywhere (with = 0), as 
well as the Taub-NUT metric as the next non-trivial example with an additional axial 
symmetry a{t) = b{t) and a removable nut singularity at the origir^^ The Atiyah- 
Hitchin metric provides an even more complicated solution, which is fully anisotropic 
and complete [1411 142| . It exhibits a removable bolt singularity at the origin, where 
b = —c and it comes asymptotically close to Taub-NUT metric with a = b as t ^ 00; 
we skip the details as they are not important for the present work. It turns out that 
these are the only complete gravitational instantons with SU{2) isometry that satisfy 
the Ricci flow equations ( 3.11| >; there is another complete metric with SU{2) isometry, 
the Eguchi-Hanson instanton, which has a = b everywhere and a removable bolt 
singularity at the origin, but its coefficients satisfy a different system of first-order 
equations. Finally, we note for completeness, that exactly the same reasoning applies 
to gravitational instantons of Einstein gravity with cosmological constant [4711331 that 
have removable nut and bolt singularities (e.g., CP^ as a gravitational instanton). 

Thus, on the one hand, in Euclidean Einstein gravity the instanton spaces have 
non-vanishing Euler number and, in many cases, they also have non-vanishing sig- 
nature (given roughly, but without any further explanation here, by the number of 



^^The power-law behavior of solutions as the flow lines approach the origin can also be read off from 
the exact solutions we have presented in various cases. 

^^This is also implied by the stability analysis around the fixed points, which shows that the time 
dependence of small fluctuations varies exponentially as t goes to infinity. 

^"^The gravitational field equations determine the Taub-NUT metric in the form shown in section 3.2, 

■-arcsinh^^ = ^^=log|^-2;« 



m mm 2m — c{t) \2m + c{t) 2m — c{t 

Setting t(, =0 for convenience, the power-law behavior of the coefficients close to the origin t = turns 
out to be fl^(t) = b'^{t) = = and describes a removable nut singularity. 
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nuts minus the number of anti-nuts | |47| that may be present). On the other hand, 
the instanton solutions of Horava-Lifshitz gravity are globally R x S'^ having zero 
Euler number and signature. They simply describe the evolution of a three-sphere 
from t = — ootof = +00 which deforms geometrically by the flow without ever be- 
coming singular along the way; as such they resemble closer the behavior of ordinary 
tnstantons in particle theories rather than the instantons of Einstein gravity. This is 
not surprising in retrospect, since consistency of the Hofava-Lifshitz gravity is not 
questionable in the projectable case for space-times with global cross-product folia- 
tion structure. In either case, the corresponding metrics are regular everywhere and 
complete and their Euclidean gravitational action is finite - though the reasoning is 
different for each theory. The finiteness of the action, which will be discussed more 
extensively shortly, makes these solutions mostly relevant in the quantum theory us- 
ing, for instance, the path integral approach. 

Let us also discuss the asymptotic structure of the solutions and compare them to 
those of ordinary gravity, since there are also important differences between the two 
theories. In Einstein gravity, the physical boundary conditions are largely determined 
by the positive-action conjecture that requires that the action of any asymptotically 
Euclidean four-metric be positive, vanishing if and only if the space is flat Il48ll33ll34| . 



Then, using the Bianchi IX form of the metric (6.2), the following possibilities arise at 



infinite proper distance t: either there is a Euclidean infinity 

a^(t) = b^{t) = c^{t) = ^t^ as t ^ ±00 (6.7) 

when < i/' < 471 (it is a conical infinity when < i/' < 2 tt) or a Taubian infinity 

a^(t) = h^{t) = t^ , = finite as i ^ ±00 (6.8) 

that encompasses the Taub-NUT metric. Combining all distinct boundary conditions 
that are available at i = and i = ±00, one ends up with a few viable solutions that 
provide the list of all complete gravitational instanton metrics with SU{2) isometry. 
Similar considerations may apply to solutions with cosmological constant. 

In Hofava-Lifshitz theory the situation is different. The Euclidean action is al- 
ways positive-definite (at least for A < 1 /3 that we are considering here) and van- 
ishes when the three-dimensional metrics are vacua of topologically massive gravity 
without any time dependence. Thus, there are no a priori conditions on the asymp- 
totic structure of instantons other than the mere existence of multiple vacua in three 
dimensions that serve as fixed points of the flow. Then, the asymptotic structure of 
space-time as ^ — > ±00 is simply determined by the specific form of the metric co- 
efficients 7i, 72 and 73 at the initial and final fixed points, respectively. Their time 
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dependence is exponential and it is completely determined by the eigenvalues of the 
characteristic matrix of small fluctuations around these fixed points. The departure 
from usual asymptotics (with zero or positive cosmological constant) is inherited to 
the solutions from the detailed balance condition and seems to be rather universal. 
It inflicts other classes of solutions, such as the construction of black-hole solutions 
whose right asymptotic structure requires departure from detailed balance using more 
general couplings |49ll50ll5T| (otherwise there is no match with observations at large 
distances). This is also closely related to the problem of obtaining ordinary gravity by 
arguing (naively) that all higher-order curvature terms are suppressed in the infrared 
regime of the theory However, it is not necessarily a big problem in the ultra-violet 
regime relevant to early time cosmology, where our discussion is applicable keeping 
A < 1/3. Abandoning detailed balance will ruin our general construction of instanton 
solutions. 



6.2 The action and moduli of SU(2) instanton metrics 

Let us now give some examples of instantons, based on the results described in previ- 
ous sections, and compute their action Sinstanton = | APV| / 2 in each case separately (see 



(2.35 1). In general, the superpotential consists of two terms PV = W^s + Weh/ which 



are given by the following expressions for Bianchi IX model geometries, 

^ (6.9) 



Wcs = 



and 



1 

1 + (^^ + 72 - 73) (71 - 72 + 73) (71 - 72 - 73) 

2717273 



1671- 



2 



Weh = 2- 



1 



V 7i7273 ^^^ + 72 + 73 - 27172 - 27273 - 27371) + 4Aw\/7i7273 

(6.10) 

and will be used next to evaluate the instanton action. 

The cases below refer to instantons constructed from interpolating trajectories of 
the Ricci-Cotton flow for different couplings, changing as the complexity of the equa- 
tions increases. 

^^For the same reason we cannot obtain the instantons of Einstein gravity from those of Hofava- 
Lifshitz theory. The first arise by dropping all higher-curvature terms and setting A = 1, whereas the 
latter exist in the full theory only for A < 1/3. Bianchi type IX models may offer a glimpse at this prob- 
lem since the gravitational potential is derived from a superpotential in both cases (see, for instance, 
p4l for the derivation of the superpotential that governs SU{2) instantons in ordinary gravity; this ref- 
erence also provides a neat qualitative picture for the completeness of their metrics using Hamiltonian 
methods). It should be easier to explore the renormalization of the coefficients of the superpotential and 
the parameter A for this class of mini-superspace models, as the theory is taken from the ultra-violet to 
the infrared domain. 
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Cotton flow. For the pure Cotton flow treated in section 3.3 there are two flxed 
points: the isotropic point with 7i = 72 = 73 = and the anisotropic fixed point 
with 7i = 72 = and 73 = L^/4x^ that arises for x = 00; two more anisotropic 

fixed points are obtained from it by permuting the three principal axes of S^. The cor- 
responding instanton, which is given in closed form by ( 3.28[ >, describes the evolution 



of a fully squashed (flattened) sphere towards the round sphere as t varies from —00 
to +00; the anti-instanton follows by reversing the time direction. 

Note that a natural entropy function associated with the volume-preserving de- 



formation of (other than W) can be defined in this case It is important for the 
mathematics of the Cotton flow, but, unlike W that determines the action of the in- 
stanton, this entropy has no deeper meaning in space-time (as far as we can tell now). 
Also, there is no (yet) known analogue of it for the combined Ricci-Cotton flow. 

The instanton solution has enhanced isometry SU{2) x 1/(1), since the deforma- 
tion line possesses axial symmetry. Also, since the interpolating trajectory is unique, 
the instanton has no moduli other than the radius of the sphere at the fixed point. In 
this case, W = Wcs and one finds that the superpotential takes the following values 
at the two fixed points, 

^iso ^ ^ ^ ^aniso ^ _ (^ -y-y^ 



Therefore, the action is 

S instanton 



47r^ 



, , (6.12) 

I Wcs I 

and it is independent of the modulus L. Obviously, there are no other instantons 
derived from the Cotton flow equations. 



^^For Bianchi IX model geometries one can define in general an additional fxinction (other than W) 
that changes monotonically under the Cotton flow. We consider 

111 

7f 72 73 

which is bounded from below by 3/(717273)^''^ = 3 ( 16 tt^/V)*''^ for a three-sphere with volume V. The 
lower bound is attained in the fully isotropic case 71 = 72 = 73. F(f ) becomes infinite when the sphere 
is completely squashed in one or more directions; as such, it is a measure of the "shape entropy" of S^. 
Using the Cotton flow | |3.18^ , we obtain 



dF 

df zfcs(7i7273)^''^ 



(72 + T3)(72-T3) , (73 + Ti)(73-7i) , (7i +T2)(7i - 72 
2 2 2 

Ti 72 73 



and, therefore, F(t) changes monotonically. For wqc, > 0, these properties oi F{t) suffice to prove the 
convergence of the flow lines to the fully isotropic fixed point regardless of initial conditions ] 
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Normalized Ricci-Cotton flow. The normalized Ricci-Cotton flow has more than 
one fixed points when ^ < 0, in which case there are instantons with SU(2) isometry 
interpolating between them as f varies for — oo to +00. 

Let us first consider the instantons connecting the two axially symmetric fixed 
points. Recall that the isotropic fixed point appears at Ji — — 73 — and 
the anisotropic point at Ji — J2 — L^/4:a and 73 = (up to permutations of the 

axes of S^), setting for notational convenience 




(6.13) 



Then, explicit calculation shows that W — Wqs + ^eh (with }i — wcsL/Ky^) takes the 
following form at the two fixed points, 

pyiso ^ ^ ^ _ 9^2 ^ Sa^AwL^) (6.14) 
wcs ^ ^ 

and ^ 

pyaniso ^ ^ (^2 + 4«^ - 14«3 + So^AwLA . (6.15) 

Wcs ^ ^ 

Although Aw does not appear in the normalized Ricci-Cotton flow equations, it enters 
into W by contributing the same at all points (recall that the volume V is preserved in 
this case). Consequently, the instanton action takes the value 

Sinstanton - ] r (« " 1)M + + 12a^ +2a + l) (6.16) 

\wcs\ ^ ' 

and it is independent of L, as expected. Notice that it vanishes when a = 1 {ji = —3), 
as required, since the two fixed point coalesce and there is no instanton in this case. 
The action is non-zero and positive for all other values ^ < 0. 

The axisymmetric solutions of the normalized Ricci-Cotton flow have been con- 
structed explicitly in section 4.4, but one should only use those branches that inter- 
polate between the two fixed points. The stability analysis performed in section 4.2 
shows that for — 3 < ^ < the isotropic fixed point is absolutely unstable and the 
anisotropic is a saddle point. Therefore, there can be only one flow line interpolat- 
ing between the two fixed points (corresponding to the axisymmetric solution we ob- 
tained) and the instanton has no moduli other than L. Exactly the same conclusion 
holds for —6 \/2 < }i < —3, since the isotropic fixed point is now absolutely stable 
and the anisotropic is a saddle point. The absence of moduli in these cases can also 
be seen schematically in figures 4 and 5, respectively. Thus, for all— 6v^2 < }i < 
the instantons have enhanced SU{2) x U{1) isometry. The situation changes dras- 
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tically when ]i < —6 \J7., since the isotropic fixed point is absolutely stable and the 
anisotropic is absolutely unstable. In this case, we have several flow lines interpolat- 
ing between the two fixed points, as can also be seen schematically in figure 6, and the 
instantons have an additional (real) modulus that labels these trajectories. The physi- 
cal interpretation of this modulus is nothing else but the geometric shape of (there 
is only one shape modulus, since the volume of space is held fixed by specifying L). 
Of course, the instanton action is independent of all moduli. 

The axisyiranetric solution we have obtained in this case has enhanced Sii(2) x 
LI(1) isometry whereas the other ones should only have an SLr(2) isometry group. 
They all correspond to regular and complete metrics on R x S^, but we have not been 
able to find them in closed form. They should be the analogue of the Atiyah-Hitchin 
metric for Horava-Lifshitz gravity when A = — oo. Their explicit construction is an 
interesting open mathematical problem. 

Finally, we turn to instantons that owe their existence to the presence of totally 
anisotropic fixed points in the problem when ^ < — 6 y/l.. They have no moduli (other 
than their volume) since they connect a saddle point with a stable or an unstable fixed 
point. These instantons also have Sli(2) isometry but no higher symmetry. 

The value of the superpotential for the totally anisotropic fixed points (see also 
section 4.2) turns out to be 

pytotal aniso = ^ ^ q + ?>a^Ky^L^) . (6.17) 

Therefore, the instanton that interpolates between these points and the totally isotropic 
fixed point has action 

72/1^ 

Sinstanton = 1 r(« + 1) • (6.18) 

Fes I 

It never becomes zero because these points cease to exist before they have the chance 
to meet with the isotropic point. Similarly, the instanton that interpolates between the 
totally anisotropic and the axially symmetric anisotropic fixed points has action 

1 C^jP- 

Sinstanton = 1 " 4) (2a3 + 1) . (6.19) 

Fes I 

The latter vanishes when = 4 (^ = — 6 \/2), as the end-points coalesce in this case, 
and it is positive definite otherwise. 

General Ricci-Cotton flow. The Ricci-Cotton flow with general couplings (provided 
that A < 1/3 and Aw is non-negative) was found to exhibit two fixed points when 
^ < 0, in which case there are instanton solutions in Horava-Lifshitz gravity. Recall 
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that the isotropic point appears at 71 = 72 = 73 = 1 / 4A w and the anisotropic point 
at 71 = 36//^/ (f/^ + 27Aw)^ and 72 = 73 = 9/ + 27Aw), assuming the presence 
of an axial symmetry 72 = 73 for all time. 

Taking j4 < 0, we define, for notational convenience, the non-negative number 

^ (6.20) 



and evaluate the superpotential W = Wcs + Weh at the two fixed points. Using ji = 
Wcs/i^w> as defined in section 5, we obtain the following results 

W^-=^(l-6.), w-=^^-^^(^--^-^^--^ + ^). (6.21) 

Therefore, the instanton action turns out to be 

12n^ (a-1)^ 



'instanton 



I Wcs I + 



(^2fl3 + 7fl2 + 24fl + 3) . (6.22) 



Note that the action is manifestly positive-definite, as required, and vanishes when 
a = 1 {}i = —3^/Aw), in which case the two fixed points coalesce and there is no 
instanton. 

Even in the presence of axial symmetry, which was used to simplify the analysis of 
the general Ricci-Cotton flow equations, we have not been able to obtain the interpo- 
lating solutions in closed form. Nevertheless, it is clear that a unique solution exists in 
this case, for all ji < 0, which interpolates between the two fixed points. For < a < 1 
{—3y/ Aw < }i < 0) the isotropic fixed point is saddle and the anisotropic is abso- 
lutely stable, and, therefore, there is a single flow line that connects the two. For a > 1 
{}i < — S^/Aw) the isotropic fixed point is now absolutely stable and the anisotropic 
is a saddle point and, therefore, the interpolating flow line is again unique. This can 
also be seen by inspecting figures 12 and 13. 

The solutions at hand have no moduli at al^^ They correspond to instantons 
with SU{2) X U{1) isometry. It will be very interesting to construct them explic- 
itly. Also, other more general solutions with strict SU{2) isometry are expected to 
exist in the general case, with 7i 7^ 72 7^ 73/ but their investigation will not be 
pursued in the present work. We only note here that all anisotropic fixed points of 
the Ricci-Cotton flow equations with general couplings seem to be axially symmetric 
even when 7i 7^ 72 7^ 73 at generic points. Thus, we expect to have instanton so- 



■^'^Even the size of at the isotropic fixed point is not free, as L was free to vary in the normalized 
flow, but it is determined by the parameters of the differential equations. Thus, it is not surprising that 
the instanton action depends on Aw (through a) in this case. 
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lutions without axial symmetry that interpolate between these fixed points, serving 
as the Horava-Lifshitz analogue of the Atiyah-Hitchin metric. They should depend 
only on one free parameter. 

The special case Aw = 0. Finally, note that as Aw is taken to zero, while keeping }i 
fixed in the general system of Ricci-Cotton flow equations, W^^° blows up to infinity, 
whereas W^^^° remains finite, tending to the value SOtz^/wqc,- Consequently, Sjnstanton 
becomes infinite and one may consider it a problem, since instantons must have finite 
action. However, in this case, there is no contradiction, since the isotropic configura- 
tion ceases to be (strictly speaking) a fixed point when it is pushed away to infinity by 
setting Aw = and, therefore, the flow line that interpolates between the two points 
(see figure 15) does not qualify as instanton solution of the theory. In conclusion, there 
are no instanton solutions when Aw = 0. 

This completes our analysis of SU{2) gravitational instantons of Horava-Lifshitz 
theory with anisotropy scaling parameter z = 3. We have obtained complete classifi- 
cation of all explicit and implicit solutions that exist for all different couplings of the 
theory satisfying the detailed balance condition, provided that A < 1 /3 and Aw > 0. 
By the same token, the variant of the theory with scaling parameter z = 2 does not 
exhibit any such instanton solutions. The results are on par with the classification of 
instantons with SU{2) isometry in Einstein gravity. The only missing technical part is 
the explicit construction of some of these instanton metrics. 

We end this section with some general remarks concerning the existence and de- 
scription of instanton metrics in Horava-Lifshitz gravity without relying on isome- 
try groups, such as SU{2). According to definition, they should be trajectories of 
the Ricci-Cotton flow equations interpolating between any two solutions of three- 
dimensional topologically massive gravity that provide the fixed points. The land- 
scape of vacua of topologically massive gravity is not known comple 
fore, it is difficult to make explicit general constructions. Also, it is rather difficult 
to investigate the general behavior of the flow equations by standard mathematical 
techniques, since they are third-order in space derivatives and even the short-time 
existence of solutions is difficult to establish in all generality. The formation and char- 
acterization of singularities is another related general open problem for these flow 
equations. Addressing these issues successfully will lead to further advances. 

^^Note, however, the recent work |52| that develops techniques to solve the field equations of topo- 
logically massive gravity (and other massive-gravity models) for three-dimensional geometries admit- 
ting a Killing vector field. Older results in this direction are neatly summarized in 1.53,1 although most 
of them focus on vacua with negative cosmological constant. 



^and, there- 
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7 Generalization to higher dimensions 



In this section we make a few remarks concerning higher-dimensional generalizations 
of Horava-Lifshitz theory and the correspondence of its instanton solutions to the 
theory of higher-order geometric flows. 



7.1 Hofava-Lifshitz gravity in 4 + 1 dimensions 

The general aspects have been reviewed in section 2 for all space-time dimensions. 
The theory is power-counting renormalizable when z = D using the appropriate su- 
perpotential W. Let us concentrate on D = 4 for definiteness, so that W is the action 
of four-dimensional gravity with higher-order corrections of the general form Q 

W[g] = I d^xy^ {aQjk,C''' + + j(R - 2Aw)) . (7.1) 

Here, Cij]^£ is the Weyl tensor and R is the Ricci scalar curvature of a four-dimensional 
Riemannian metric g that describes the geometry of spatial slices in a five-dimensional 
space-time with topology Ais = R x Ad/^. Here, there is no need to include the term 
RijR^^ because it can be removed by a Gauss-Bonnet topological term, adjusting the 
coefficients a and fi. 

Thus, Horava-Lifshitz gravity in 4 + 1 dimensions with anisotropic scaling z = 4 
is defined by the action 

(7.2) 

using the extrinsic curvature Kij of Ai^ and the metric G^j^^ of superspace with param- 
eter A. Also, following the general discussion of section 2, we will also take A < 1 /4 
so that the Euclidean counterpart of this action is manifestly bounded from below. 

The theory with detailed balance is completely specified by the choice of W[g] . It is 
given by the general expression ( 7.1[ > in D = 4; other appropriate choices of W should 



be made in higher dimensions to render the theory power-counting renormalizable. 
We also note for completeness that if higher-order curvature functionals are chosen in 
D spatial dimensions so that z > D, the resulting gravitational theory will be power- 
counting superenormalizable [ZJ- Such generalizations will not be considered at all in 
the present work. 

Next, we illustrate the structure of the resulting equations by considering the sim- 
plest higher-dimensional case with z = D = 4. 
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7.2 Bach flow and its variants 



Solutions of the Euclidean five-dimensional Horava-Lifshitz gravity can be obtained 
from the geometric-flow equation 

^8*, = ±^V7j^ + v,e,+v,e, (7.3) 

that describes deformations of the four-dimensional Riemannian metric gij. The lapse 
function N{t) can be set equal to 1 by time redefinition. 

The details can be worked out using the following identity, which is well known 
among people working in conformal Weyl gravity, 

Weyl^ (7.4) 

where 

Wweyl = / d'x^QjkiC'j^' (7.5) 
is the quadratic Weyl tensor action functional and 

B^j = Vk^iC^^' + IrmC''' (7.6) 

is the so called Bach tensor ||54| . It is a fourth-order symmetric and traceless tensor 
that clearly vanishes when the four-dimensional metric is conformally flat. The Bach 
tensor provides the analogue of the Einstein tensor in the field equations of conformal 
Weyl gravity, and, as such, it is also covariantly conserved. 

Thus, for this particular choice of superpotential W, the corresponding geometric 
flow takes the form ^ 

Mrj = T\B,j + V,^j + Vj^,, (7.7) 

and it can be naturally called Bach flow. Its fixed points (modulo reparametrizations) 
are the vacuum solutions of conformal Weyl gravity and include the isotropic (con- 
stant curvature) metric on S^. It is mathematically more interesting to pick the sign 
that drives the evolution towards the fixed points rather that away from them. Al- 
though this is a higher-order flow, it is better behaved mathematically than the third- 
order Cotton flow. Thus, one should investigate it in detail and attempt to construct 
solutions. It is a new geometric flow that has not appeared in the mathematics liter- 
ature before, to the best of our knowledge. By restricting it to Kahler manifolds, it 
might be also interesting to compare it (and the variants which are discussed below) 
with other well known geometric flows of fourth-order, such as the Calabi flow ||55| . 
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If there is an additional contribution to W given by the quadratic Ricci scalar cur- 
vature action, 

Wr2 = J dV^R2^ (7.8) 
it will account for the gradient term 



(7.9) 



with 



W = 2g'j\/kS/^R - IV'VJR - 2RR'j + ^'^R^ . (7.10) 

This tensor is symmetric but not traceless. Then, the complete flow equation will be 
a variant of the Bach flow receiving contributions from B'^ and H'^, which are both 
fourth-order. Of course, in the general case, there will also be subleading curvature 
terms associated with the Einstein tensor G'^ by adding the four-dimensional Einstein- 
Hilbert action (possibly with a cosmological constant) to the superpotential W. 

Instanton solutions will correspond to flow lines interpolating between different 
vacua of four-dimensional conf ormal Weyl gravity (and its deformations thereof), but 
again it seems rather difficult to derive explicit general results. Using four-dimensional 
model geometries may provide some simple and tractable mini-superspace models 
that are worth studying in the future. 

Similar considerations apply to all higher-dimensional generalizations of Horava- 
Lif shitz gravity. In D + 1 dimensions, the non-relativistic gravitational theory becomes 
power-counting renormalizable when the anisotropic scaling parameter is z = D. 
Then, for the appropriate choice of W, we obtain geometric flows of order D that 
describe instanton-like configurations of the Euclidean (D + l)-dimensional theory 
when A < 1/D. This framework hosts very naturally a whole hierarchy of geometric 
flows and provides a reason to study them. 



8 Conclusions and discussion 

We examined the Euclidean version of Horava-Lif shitz gravity satisfying the detailed 
balance condition and described its instanton solutions as flow lines interpolating be- 
tween different fixed points of a new class of geometric evolution equations, which 
are first-order in time. Although the specific couplings implied by detailed balance 
are rather restrictive (and sometimes appear to be problematic), the general connec- 
tion between instanton solutions and geometric flows is rather interesting in many 
respects. Focusing to 3 + 1 dimensions, where the potential term is derived from a su- 
perpotential W given by the action functional of three-dimensional topologically mas- 
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sive gravity and the anisotropy scaling parameter of the theory is z = 3, the driving 
curvature terms are provided by a certain combination of the Cotton and Ricci tensors 
as well as the cosmological constant term. The geometric-flow equations, called Ricci- 
Cotton flow, were shown to exhibit an entropy functional that is given by W and can 
be used to put a lower bound on the Euclidean Hof ava-Lifshitz gravitational action. 

Our construction requires A < 1/3 and Aw > 0, but otherwise the parameters of 
the theory can be arbitrary within the class of detailed balance couplings. Fixed points 
of the flow are provided by classical solutions of the three-dimensional topologically 
massive gravity and they correspond to static solutions of the (3 + l)-dimensional the- 
ory. As such, they include constant-curvature isotropic metrics in three dimensions 
as well as anisotropic configurations obtained by balancing the deformation effect of 
the Cotton and Ricci tensors. Since there is no general classification of these metrics, 
the landscape of fixed points remains largely unexplored to the best of our knowl- 
edge. Running solutions represent genuine time-dependent configurations, but they 
are even more difficult to investigate in exact terms. Thus, the Ricci-Cotton flow ap- 
pears to be a rather complex system of equations that deserves proper mathematical 
study on general grounds. Addressing these problems in all generality remains out of 
reach at the moment, but some simple mini-superspace truncations of the equations 
help to obtain concrete results in simple cases that are interesting both physically and 
mathematically. 

We found that the homogeneous model geometries on three-manifolds provide 
consistent truncation of the Ricci-Cotton flow equations. In particular, focusing on 
the Bianchi IX model geometries on S^, so that the corresponding gravitational instan- 
tons exhibit SU{2) group of isometries, we were able to classify the fixed points of 
the flow (isotropic as well as anisotropic) and study their stability properties for a va- 
riety of different couplings. Some special solutions with axial symmetry (associated 
with SU{2) X U{1) isometry group) were constructed explicitly and their space-time 
interpretation was discussed in analogy with the gravitational instanton solutions of 
ordinary gravity. In particular, we have arrived at complete classification of the in- 
stanton solutions with SU{2) isometry. It remains to be seen whether more general 
running solutions can be constructed explicitly beyond their qualitative description 
based on the phase portraits of the flow. Also, it will be interesting to find other con- 
sistent reductions of the flow equations beyond the class of homogeneous geometries, 
but we have not yet been able to obtain any concrete results in this direction. 

Another possibility that has not been discussed at all in this paper is to consider 
super-renormalizable versions of Hofava-Lifshitz gravity in 3 + 1 dimensions with 
anisotropic scaling z = 4. These are generated by a superpotential W - other than the 
action functional of topologically massive gravity - which contains higher-order Ricci 
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curvature terms such as RijR^^ and on top of the cosmological Einstein-Hilbert ac- 
tion in three dimensions |lZ|- In this case, the instanton solutions will be described 
by geometric flows in three dimensions with fourth-order derivatives in their driv- 
ing curvature terms. The resulting equations appear to have some nice mathematical 
properties (compared to the third-order Ricci-Cotton flow) and they also seem to ad- 
mit consistent reduction to an autonomous system of ordinary differential equations 
for homogeneous model geometries. A particular choice of such W is provided by 
the action of the so called "new massive gravity" in three dimensions that contains 
both terms RijR^^ and J^^ with relative coefficient —3/8 [|56il . We intend to investigate 
elsewhere the corresponding fourth-order flows ||57|, together with the associated in- 
stanton solutions, and examine the privileged role (if any) of new massive gravity in 
this context. 

Higher-dimensional generalizations were also briefly discussed. It was pointed 
out that instanton solutions exist in all dimensions and their defining equations pro- 
vide new classes of geometric-flow equations, such as the Bach flow in four dimen- 
sions. In general, the driving curvature terms of such flows contain spatial deriva- 
tives of order z (equal to the anisotropy scale parameter that renders the higher- 
dimensional Horava-Lifshitz gravity power-counting renormalizable) and they de- 
scribe metric deformations on spatial slices of dimension D = z. The hierarchy of 
such flows has not been considered in the literature before and they certainly pose 
several interesting questions that are worth studying in the details. They should also 
be of interest to the mathematics community working on geometric analysis. In all 
cases, the non-relativistic theory of gravity provides a general framework to embed 
geometric evolution equations. The situation should be compared to general rela- 
tivity and string theory, where such embedding is only possible in some very special 
cases, such as the Ricci flow on homogeneous three-geometries that can be interpreted 
as self-dual gravitational instantons in four dimensions or using some appropriately 
chosen higher-dimensional plane-wave gravitational backgrounds. 

The off-shell formulation of string theory based on the world-sheet renormaliza- 
tion group equations provides a natural framework for the appearance of the Ricci 
flow (and other closely related geometric -flow equations) in gravitational physics. In 
this context, closed string tachyon condensation is described by transitions from one 
fixed point to another more stable fixed point. Thus, the lines of the Ricci flow resem- 
ble instanton transitions among different vacua of the string landscape. In Hofava- 
Lifshitz gravity, on the other hand, the Ricci flow (and its variants) describe the in- 
stantons of the theory when the anisotropic scaling is z = 2. Therefore, it seems in- 
teresting to investigate further this aspect while searching for possible embedding of 
the non-relativistic theory of gravity into a more fundamental theory. Likewise, non- 
relativistic theories of gravity with higher anisotropic scaling, in particular z = D, and 
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their instanton solutions may admit a similar description and interpretation in terms 
of a more fundamental theory. It is not yet known, however, whether the geometric 
evolution equations we are considering here can also arise as renormalization-group 
equations in a class of quantum field theories. 

Finally, another interesting problem is the use of instantons for the quantization of 
Horava-Lifshitz gravity. One possible line of work in this direction is the path integral 
approach over Euclidean space-times with applications to quantum cosmology in the 
spirit of Hartle-Hawking proposal. The quantization of mini-superspace models ap- 
pears to be tractable, at least for homogeneous (but generally non-isotropic) geome- 
tries, and requires special attention. They can also provide some non-perturbative 
information about the quantum theory and a testing bed for comparison with the 
quantization of ordinary gravity. 

It remains to be seen whether the non-relativistic theory of gravity is a viable al- 
ternative to Einstein gravity at very short distances. However, the simplified version 
of the theory with detailed balance can also play another role in physics, serving as 
landscape explorer of the vacuum structure of relativistic field theories determined by 
W (with topologically massive gravity being just an example). It provides an effective 
particle model to describe transitions among different vacua through instantons. It 
also offers a dynamical principle for vacuum selection that is worth exploring further 
in all generality using the powerful tools of geometric flows and associated entropy 
functions. Advocating this point of view introduces a new twist to the subject and de- 
parts from the idea (and the problems that seem to accompany it) that Hof ava-Lif shitz 
gravity is the ultra-violet completion of a fundamental theory. It could have also been 
used from the very beginning as an alternative motivation for the present work. 

We hope to return to these topics elsewhere in the near future. 
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